Bond-order and charge-density waves in the isotropic interacting two-dimensional 
quarter-filled band and the insulating state proximate to organic superconductivity 
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We report three surprising results regarding the nature of the spatial broken symmetries in the 
two-dimensional (2D), quarter-filled band with strong electron-electron interactions that provides 
a microscopic model of the 2:1 cationic organic charge transfer solids (CTS). First, in direct con- 
tradiction to the predictions of one-electron theory, we find a coexisting "bond-order and charge- 
density wave" (BCDW) insulating ground state in the 2D rectangular lattice for all anisotropies, 
including the isotropic limit. Second, in contrast to the interacting half-filled band, which ex- 
hibits one singlet-to-antiferromagnet (AFM) transition as the interchain coupling is increased from 
zero, there occur in the interacting quarter-filled band two distinct transitions: a similar singlet-to- 
antiferromagnet/spin-density wave (AFM/SDW) transition at small interchain coupling, giving rise 
to a bond-charge-spin density wave (BCSDW) state, followed by a second AFM/SDW-to-singlet 
transition at large interchain coupling. Third, we show that our conclusions remain unchanged if 
one assumes the conventional "effective 1/2-filled" lattice of dimer sites for the CTS: the dimer lat- 
tice unconditionally dimerizes again to give the same BCDW found in the quarter-filled band. We 
make detailed comparisons to recent experiments in the tetramethyl-tetrathiafulvalene (TMTTF), 
tetramethyl-tetraselenafulvalene (TMTSF), bisethylenedithio-tetrathiafulvalene (BEDT-TTF) and 
bisethylenedithio-tetraselenafulvalene (BETS)-based CTS. Our theory explains the mixed charge- 
spin density waves observed in TMTSF and certain BEDT-TTF systems, as well as the absence of 
antiferromagnetism in the BETS-based systems. An important consequence of this work is the sug- 
gestion that organic superconductivity is related to the proximate Coulomb-induced BCDW, with 
the SDW that coexists for large anisotropies being also a consequence of the BCDW, rather than 
the driver of superconductivity. We point out that the BCDW and BCSDW states are analogous to 
four different classes of "paired" semiconductors that are obtained within certain models of exotic 
superconductivity. That all four of these models can in principle give rise to superconductivity in 
the weakly incommensurate regime provides further motivation for the notion that the BCDW may 
be driving the superconductivity in the organics. 



I. INTRODUCTION 

Theoretical discussions of spatial broken symmetries 
in strongly correlated electron systems have largely fo- 
cused on the 1 /2-filled band Mott- Hubbard semiconduc- 
tor. The one-dimensional (ID) case has n been widely 
discussed in the context of polyacetylenea-B. Here it 
is known that Coulomb electron-electron (e-e) interac- 
tions can strongly enhance the 2kp (kp = one-electron 
Fermi wavevector) bond-alternation expected within the 
Peierls purely electron-phonon (e-ph) coupled model, giv- 
ing rise to a periodic modulation of the bond-order, a 
bond-order wave (BOW). In the limit of very strong on- 
site Coulomb interaction, the BOW instability is usu- 
ally referred to as the spin-Peierls (SP) instability. In 
the presence of intersite Coulomb interactions, and for 
certain relative values of the on-site and intersite inter- 
action parameters, a charge-density wave (CDW), peri- 
odic modulation of the site charge density, can be the 
dominant instability 3 ]. The BOW and the CDW occur 
in largely nonoverlapping regions of the parameter space 



and compete against each otheioQ. True antiferromag- 
netism (AFM) — ie, a long-range order (LRO) 2kp spin- 
density wave (SDW) — is absent in for spin-rotationally 
invariant models in ID, and the ground state is domi- 
nated by singlet spin coupling, which favors the BOW 
over the SDW. Two-dimensionality is thus essential for 
the SDW. 

The 1/2-filled isotropic two-dimensional (2D) case has 
been investigated in great detail in recent years (mostly 
for the case of large inttasite Hubbard interaction but 
zero intersite interaction)?!, as this limiting case is known 
to describe the parent semiconductor compounds of 
copper-oxide based high temperature superconductors. 
The BOW instability that characterises the ID chain is 
destabilized in 2D by Coulomb interactiorffil3, and the 
dominant broken symmetry here is the 2hp SDW, with 
periodic modulation of the spin density. Most recently, it 
has been demonstrated that this SDW state appears for 
the smallest nonzero interchain hopping in weakly cou- 
pled 1/2-filled band chainsEj, in agceeaaent with previous 
renormalization group calculationstSO. As in 1DD, there 
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is no CDW-SDW coexistence in The absence of 

coexistence between the BOW and SDW for the 1/2- 
filled band in both ID and 2D can be readily understood 
intuitively: the BOW requires spin-singlet coupling be- 
tween alternate nearest neighbor spins, which clearly has 
to disappear in the SDW. An alternate way of viewing 
this is to observe that the probability of charge-transfer 
to the left and to the right in the AFM are exactly equal, 
and therefore the SDW cannot coexist with the BOW. 
On the other hand, both the BOW and the SDW re- 
quire that the site-occupancies by electrons are strictly 
uniform, and thus neither the ID BOW nor the 2D SDW 
will coexist with the CDW. 

Coupled 1/2-filled band chains have also been dis- 
cussed r-within the context of the so-called ladder 
systemsE-i Whether or not a given n-leg ladder system, 
for small n, exhibits the BOW now depends on whether 
n is odd or even. This feature of the ladder systems could 
have been anticipated from the physics of the odd versus 
even S Heisenberg chainstj. Thus at least for the simplest 
monatomic lattices, ground states of the 1/2-filled band 
are known: the BOW, CDW and SDW phases compete 
against one another and do not coexist, and 2D behavior 
emerges for the smallest 2D coupling. 

In contrast to the 1/2-filled band, broken symmetries 
in non- 1/2-filled bands with strong e-e-interactions have 
been investigated primarily in ID limitllTcj or at_most in 
the quasi-lD regime of weak interchain couplingEa. This 
emphasis likely arises from the theoretical preconception 
that finite one-electron hopping between chains destroys 
the nesting feature that characterizes the ID limit, lead- 
ing necessarily to the restoration of the metallic phaseO. 
A recent work has examined coupled chains in the limit 
of weak e-e interactionsE^I. The weak-coupling approx- 
imation employed in reference [23] reproduces the loss 
of nesting predicted within band theory. While the con- 
tinuum renormalization group calculationalj'EJ predicted 
CDW-SDW coexistence for incommensurate bandfillings, 
early quantum Monte Carlo calculations for the 1 /4-fillcd 
band failed to find this coexistenceL3. Many more re- 
cent numerical simulations on discrete finite systems as- 
sume the absence of coexistence between the 2kp BOW, 
the 2kp CDW and the 2kp SDW that characterizes that 
1/2-filled band also applies to the non- 1/2-filled bands. 
Indeed, it is often assumed that the CDW is driven by 
the e-ph interactions and the SDW by e-e interactions 
and that their effects are competing. This assumption is 
made despite the result mentioned above that already in 
the simplest case of the ID 1/2-filled band, e-e and e-ph 
interaction effects are known not to be competing but 
to act. in a co-operative way to give the enhanced 2kj? 
BOWtl. 

Recently, we have begun a systematic study of the na- 
ture of the broken symmetry ground states in the 2D 
1 / 4-filled band on an anisotropic-tactangular lattice with 
both e-ph and e-e interactionaHJo. Earlier work by us 
had already established the cooperative coexistence be- 
tween the BOW and the period 4 "2k F " CDW in the 
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FIG. 1. The "effective" half-filled band and dimerized 
dimer model, (a) Dimerization in a 3/4-filled ID band of 
electrons leads to a gap in the single particle spectrum at 
k = ±7r/2a (a = lattice spacing), resulting in a half-filled 
upper subband. Note that although the actual CTS mate- 
rials are indeed nominally 3/4-filled electron bands (hence 
1/4- filled hole bands), in the text we follow the convention 
and refer to them simply as "1/4- filled", (b) A real space 
depiction of a 2D lattice of dimers in the strong correlation 
limit. The two sites within the parentheses form one lat- 
tice point of the dimer lattice, and the intradimer bonds are 
stronger than the interdimer bonds. The charge and spin 
populations on individual sites within each dimer are equal, 
and the effective 1/2-filled band lattice is antiferromagnetic 
in 2D. (c) Schematic of a frozen valence bond state resulting 
from the dimerization of dimer lattice. The interdimer bonds 
are now different; the line denotes a singlet bond. This frozen 
valence bond diagram is relevant in the ID limit, and then 
again for the strongly 2D case, where the antiferromagnetism 
has disappeared. The antiferromagnetic phase that occurs for 
intermediate interchain coupling is shown in Fig. 2. 

ID 1/4-filled band, with each broken syipsnetry 
enhancing r-thc other, for both nonintcractingc3 and 
interactingEO electrons. The latter results have been sub- 
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sequently confirmed- by Riera and PoilblanccS In the 
more recent workcju3 we have demonstrated an appar- 
ently unique feature of the 1/4-filled band: namely, the 
coexistence of the BOW-CDW with the period 4 "2k F " 
SDW, giving rise to a coupled Bond-Charge-Spin density 
wave (BCSDW) that appears for weak interchain electron 
transfer between chains. 

In the present paper, we extend our calculations to the 
full range of anisotropics, from uncoupled chains to an 
isotropic 2D lattice. We include both the SSH intersite 
phonons that driv|8-,a BOWtl and the Holstein phonons 
that drive a CDWcp. We list three primary motivations 
for this extensioncJ. First, the cooperative coexistence 
between the BOW and the 2k F SDW found in the 1/4- 
filled band for weak interchain transfer is exactly opposite 
to the competition between the 2kp BOW and the 2kp 
SDW (with the latter dominating for nonzero interchain 
transfer) in the 1/2-filled band. It is then immediately 
natural to ask what the nature of the ground state is for 
strong interchain hopping of electrons in the 1/4- filled 
band. Second, from a more general theoretical perspec- 
tive, whether or not the vanishing of density waves that is 
predicted by one-electron nesting ideas remains true for 
strongly correlated electrons is of considerable general in- 
terest. Finally, our results are likely to have relevance to 
experimental observations in the orgMic-CTS, including 
those that exhibit superconductivityE2l~E2l. 

Our investigations yield the surprising result that the 
coexisting Bond-Charge density wave (BCDW) persists 
as the ground state of the strongly correlated 1/4- filled 
band in 2D for all values of the interchain electron trans- 
fer, including the isotropic limit. We show that this 
result can be understood physically as a consequence 
of interchain confineraeni,arising from strong intrachain 
Coulomb interactionsE3~t3. The SDW component of the 
BCSDW, on the other hand, attains a maximum am- 
plitude at some intermediate interchain transfer, after 
which it typically vanishes at a critical value of the trans- 
fer. 

In order to discuss applications of results to real mate- 
rials, including the 2:1 cationic CTS, we need to clarify 
an important aspect of our approach vis-a-vis most pre- 
vious work on models of these materials. In our above 
discussion of band-filling, "1/4-filled" is defined in the 
usual manner: namely, in the absence of the BCDW, the 
lattice is uniform in at least one direction, and the aver- 
age density of electrons per site is 1/2. In real materials, 
crystal structure effects often cause a lattice dimerization 
that is unrelated to any underlying electronic or mag- 
netic instability (see below)E3. As shown in Fig. 0(a), 
this dimerization leads to a gap in the single electron 
spectrum at fcp = n/2a, and consequently suggests us- 
ing an effective 1/2-filled band model that focuses on the 
upper subband. In real space terms, this approximation 
amounts to considering the system as a set of (tightly 
bound) dimers (i.e., a diatomic lattice) with one electron 
per dimer site, as showH-jn-Jigure. [l](b). This approach 
has been widely appliedE3c3, particularly with consider- 



able success in the context of the magnetic field-induced 
spin density wave (FISDW) in 2:1 salts of TMTSFBEI. 
As we show below, a further dimerization of the dimer 
lattice is unconditional in both ID (the well-known spin- 
Peierls transition) and 2D (a surprising new result), and 
that this dimerization of the dimer lattice leads spon- 
taneously to different electronic populations on the sites 
within a dimer, i.e., to the same 2Uf CDW that occurs 
in the 1/4-filled (monoatomic) band (see Fig. 1(c)). For 
small interchain electron transfer, the BCSDW will there- 
fore have nearly the same structure as the original 1/4- 
filled band. This is a third new result, perhaps also sur- 
prising, and shows that the number of electrons per site 
within a unit cell is a more fundamental parameter than 
the bandfilling: the latter is strictly a one-electron con- 
cept of limited use in the interacting electron picture. 

We expect our results to be relevant for the ID semi- 
conductors (TMTTF) 2 X, the so-called "quasi-lD" or- 
ganic superconductors (TMTSF^X, as well as the 2D 
organic superconductors (BF,D T=T TF)nX and the more 
recently synthesized (BETS^Xcil In reference [24] we 
showed-jpthiat the highly unusual "mixed CDW-SDW 
state"Ny>B found in (TMTTF) 2 Br, (TMTSF) 2 PF 6 
and (TMTSF)2AsF6 can be explained naturally as the 
BCSDW state within the strongly correlated 1/4-filled 
band scenario. Our current work shows that dimeriza- 
tion of the dimer lattice leads to the same results, and 
hence therweak high temperature dimerization along the 
stack axisE3 is effectively irrelevant: starting from either 
the 1/4-filled model or the effective 1/2-filled scenario, 
the final outcome is the samec3. 

With these comments complete, we can describe the 
organization of the remainder of the paper. In Section 
II we introduce our model Hamiltonian, as well as that 
of the dimerized dimer model. In Section III, we present 
physical, intuitive arguments, based |-tp.[a_ configuration 
space picture of broken symmctryoQtSEH that predict 
both the BCSDW for weak interchain electron transfer 
and the persistent BCDW state in the isotropic limit. In 
Section IV, we present the results of extensive numerical 
studies, exploring behavior in both the strict ID limit 
and for the full range of anisotropics in the quasi-2D 
case. These studies, in confirmation of the qualitative 
predictions of Section III: (i) establish the persistence of 
the BCDW up to the isotropic limit; (ii) suggest the oc- 
currence of two quantum critical transition as an SDW 
first appears for weak transverse hopping and then dis- 
appears for the nearly isotropic case; and (iii) prove the 
equivalence of the 1/2-filled dimerized dimer and 1/4- 
filled monatomic lattices. For clarity, in Section V we 
summarize our theoretical conclusions; readers not inter- 
ested in the underlying physical arguments or numerical 
details can skip directly to this summary in Section V. 
In Section VI, we examine in some detail several recent 
experiments that indicate the applicability of our theory 
to the insulating states that are observed to be proxi- 
mate to the superconducting states in the organic CTS. 
Finally, in Section VII, we indicate possible future di- 
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rections for our research, focusing on commensurability 
defects in the BCDW state and their possible role in the 
proximate superconducting phases. We point out sev- 
eral intriguing similarities between this potential micro- 
scopic mechanism for superconductivity and other recent 
phenomenological models. We conclude the article with 
three appendices, which deal with various more technical 
arguments and details of the numerical methods. 

II. MODELS AND OBSERVABLES 

We consider two different extended Peierls-Hubbard 
Hamiltonians on a rectangular lattice 2D with (in gen- 
eral) anisotropic electron hopping. The first model de- 
scribes a monatomic 1 / 4- filled band and is defined by the 
Hamiltonian 

H = Hq + H ee + Hi n t er (la) 

Hp = - [t - a(Aj tM )]Bjj + i <M ,M,a + J~] Vj,M n J,M 

j,M,a j.AI 

+K t /2 £(A- M f + K 2 /2 J2 v j,M (16) 

j,M j,M 

H ee = U^2 n j,M^n jtM ,l + V 2J Hj,M«j+l,M (1c) 
j,M j,M 

Hinter = ~tx ^ ^i,j,M,M+l,<T (1^) 

In the above, j is a site index, M is a chaip-iiftdex, a is 
spin, and we assume a rectangular lattice^OcZl. As ij_ 
varies from to t, the electronic properties vary from ID 
to 2D. An implicit parameter in the above Hamiltonian 
is the bandfilling, or more precisely p. We shall focus on 
the 1/4-filled case, for which p = 1/2. In applications to 
the organic CTS, each site is occupied by a single organic 
molecule, the displacement of which from equilibrium is 
described by u jtM (with A jtM = (u j+ i tM ~ Wj.JVf)); Vj,M 
is an intra-molecular vibration, rij,M,o — cj M a Cj t M,a, 

n i,M = J2a n 3,M^> and Bj,k,L,M,<r = i c ].L.a C k,M ,o~ + h.c] , 

where cj L a is a Fermion operator. We treat the phonons 
in the adiabatic approximation and are interested in un- 
conditional broken symmetry solutions that occur for e- 
ph couplings (a,/3) — > + . All energies such as U, V, 
and t± will be given in units of the undistorted intra- 
chain hopping integral t. 

The second model describes a diatomic/ dirtier lattice, 
with one electron per dimer. The Hamiltonian for this 
case is similar to that above, with identical H ee and 
Hinter, but with modified intrachain one-electron term 



Hq = —tl 2_J B 2 j-l,2j,M,ALa (2) 
3,M,<r 

— ^ \p2 ~ OtAj,Af]B2j,2j+l,M,M,<7 + ~W ^ (^j,M Y 
j.M,a j,M 

In the above each pair of sites (2j-l,M) and (2j,M) 
forms a dimer with fixed hopping t± > t between them, 
Aj )Af = {u 2 j+\M ~U2j,m), with U2j-i,M = U2j x M\ this 
means that there is no modulation of the intradimer bond 
length, and the dimer unit is displaced as a whole. As 
written, the model assumes an "in-phase" 2D arrange- 
ment of the dimer units (i.e., dimers on different chains 
lie directly above one another), which we have deter- 
mined to be the lower energy configuration for both zero 
and nonzero A^m- Notice that H' does not contain 
the Holstein on-site e-ph coupling. Nevertheless, we will 
show that a site-diagonal CDW is a consequence of the 
BOW here. 

The broken symmetries we are interested in are (i) 
the BOW, with periodic modulations of the intrachain 
nearest neighbor hand order Bj,j+i,M,M,a)] (h) the 
CDW, with periodic modulations of the site charge- 
density (nj M)i an d (hi) the SDW, with periodic mod- 
ulations of the site spin-density (rij t M,i — n j,M.i)- Note 
that in case of the dimer lattice (Eq. (2)) we are inter- 
ested in both intra- and interdimer charge and spin mod- 
ulations, although bond modulations can occur only be- 
tween dimers. Furthermore, in the CDW and the SDW 
the modulations of the site-based densities occur along 
both longitudinal and transverse directions (though not 
necessarily with the same periodicities, see below). In 
case of the BOW, a complete description would require 
the determination of the phase difference between con- 
secutive chains. 



III. CONFIGURATION SPACE PICTURE OF 
SPATIAL BROKEN SYMMETRY 

The physical arguments presented in this section pro- 
vide crucial insights that allow us to anticipate the appar- 
ently counterintuitive results of this paper. The need to 
develop such arguments arises from the limitations inher- 
ent in all true many-body numerical simulations of strong 
correlated electron systems: namely, one can study only 
systems of limited size and distinguishing finite-size ar- 
tifacts from true results requires physical understanding. 
In turn, true many-body numerical methods are essen- 
tial here because of the intermediate magnitude of the 
e-e interactions (comparable to the bandwidths) in the 
organic CTS, which renders both mean field and pertur- 
bation theoretic approaches questionable. For instance, 
even in the strictbi.lD limit, where well-established RGt£l 
and bosonizationO techniques have existed for decades, 
for the intermediate coupling regime, there have recently 
been some surprising discover ies , in the phase diagram of 
the extended Hubbard modelE§c3. In 2D, developing a 
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clear physical intuition is still more crucial, as numeri- 
cally tractable lattices are even farther from the thermo- 
dynamic limit, and the competition among broken sym- 
metries is likely to be more subtiy.-J3rief presentations of 
these physical ideas for t±_ = (J23'E3 and tj_ << £-3 have 
been made previously. Here we discuss these ideas for the 
complete range < tj_ < t, focusing on (i) the transition 
from ID to 2D, and (ii) the difference from the 1/2-filled 
band monatomic lattice. 

A physical picture of spatial broken symmetry in 
strongly correlated electron systems must necessarily 
be based on configuration space ideas, as one-electron 
bands have simply ceased to exist for strong e-e inter- 
action. WithmJic configuration space picture of bro- 
ken symmetryB'ErEj, each broken symmetry state, inde- 
pendent of band-filling, can be associated with a small 
number of equivalent configurations that are related by 
the symmetry operator in question. For commensurate 
p, these configurations are easily determined by inspec- 
tion. The relevant configurations consist of repeat units 
which themselves possess the same periodicity as the 
density wave. For illustration, we choose the ID 1/2- 
filled band. In this case, each broken symmetry has two 
extreme configurations, the pairs corresponding to the 
SDW, BOW and CDW being, respectively: the two Neel 
states ... tltl ■■■ and ... ||IT •■• (SDW); the two near- 
est neighbor valence bond diagrams (1,2)(3,4)(5,6)....(N 
- 1,N) and (N,1)(2,3)(4,5)....(N - 2, N - 1) (where (i,j) is 
a spin singlet bond between sites i and j and N is the num- 
ber of sites) (BOW); and the configurations ...202020... 
and ...020202... (where the numbers denote site occupan- 
cies) (CDW) . N applications of the one-electron hopping 
term in Eq. (I) on any one extreme configuration (cor- 
responding to a given broken symmetry) generates the 
other extreme configuration, but for N — > oo this mix- 
ing of configurations is small, and the ground state re- 
sembles one or the other of the extreme configurations 
qualitatively, with reduced spin moment, bond order _pr 
charge-density difference due to quantum fluctuations! 

The key insight of the configuration space heuristics 
is that the qualitative effects of many-body Coulomb in- 
teractions, as well as additional one-electron terms, can 
be deduced frampthcir effects on any one of the extreme 
configurationsmlEl As a trivial example of this, a repul- 
sive Hubbard U destroys the CDW in the 1 /2-filled band, 
simply because double occupancies in the extreme con- 
figuration ...202020... "cost" prohibitively high energy. 
Significantly, in the I/2-filled band, the extreme config- 
urations favoring the SDW, the BOW and the CDW are 
different, and there is a complete lack of overlap between 
them. This essentially guarantees the absence of coex- 
istence among these broken symmetries in both ID and 
2D. 

To apply these ideas to the ID 1/4-filled band, we be- 
gin by considering the on-site charge configurations. A 
2kp (4kp) density wave here has period 4 (2) in config- 
uration space. As discussed above, the extreme config- 
urations of interest must also have period 4 or 2, and 



there are then only three distinct sets of extreme charge 
configurations. These contain the repeat units ...2000..., 
...1100..., and ...1010..., respectively, where the numbers 
again denote site occupancies. There are four distinct 
configurations for sets 1 (...2000...) and 2 (....1100....), 
whereas there are only 2 for set 3 (....1010....). By anal- 
ogy with the 1/2-filled band (see above), we now intro- 
duce spins and note that configurations belonging to sets 
2 and 3 can again have spin singlet bonds between pairs 
of nearest neighbor singly occupied sites, or the spins of 
the occupied sites can alternate as in the 1/2-filled band 
Neel configurations. Let us now show, by considering 
the different cases separately, how e-e interactions affect 
these configurations and how an understanding of these 
effects suggests (correctly!) the broken symmetries to be 
studied. 

A. 1/4-filled band, t± = 0, U = V - 

The non-interacting case provides a simple example 
to introduce some of the important differences between 
the 1/4-filled and 1/2-filled bands. Actual calculation 
indicates that within the ID Holstein model the-xharge 
densities pj on the sites have the functional formE£l 

Pj = 0.5 + po cos(2kp ja) = 0.5 + po cos(nj/2) (3) 

This charge density pattern could have been anticipated 
by focusing on the extreme configuration ...2000..., which 
also predicts three different charge densities (large, inter- 
mediate, small and intermediate), since each '0' that is 
immediately next to a '2' is different from the other pair 
of sites labeled '0' that are further away from the '2'. 
occupancy scheme ...2000..., the probabilities of charge- 
transfer between a '2' and the two neighboring '0's are 
larger than that between the two neighboring '0's them- 
selves. For nonzero a in Eq. (1), this difference in charge- 
transfers leads to lattice distortion of the form 

uj = uocos(2kpja) — uq cos(7t//2), (4) 

with bonding pattern "SSWW" (for strong, strong, weak, 
weak), where a strong (weak) bond has hopping tg > t 
(tw < t) . This then is one very important difference from 
the 1/2-filled band: whereas in the 1/2-filled band dif- 
ferences in bond-orders arise from spin-effects only (the 
probability of charge-transfer is greater between near- 
est neighbor singlet-couplad sites than between nearest 
neighbor non-bonded sitestil), in non- 1/2-filled bands this 
difference can also originate from site occupancies. Pre- 
cisely because the BOW and the CDW here are both de- 
rived from the same extreme configuration, they coexist 
in the noninteracting 1 / 4- filled bandEa 

B. 1/4-filled band, t± = 0, U, V > 

For nonzero (positive) U and V, the interplay among 
the various possible broken symmetries becomes both 
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more subtle and more interesting. Since double occupan- 
cies "cost" energy, the extreme configuration ...2000... 
is suppressed even at a relatively small [/E3. For the 
strongly correlated (U — + oo) ID 1/4-filled band with 
convex long range interactions, Hubbard showed that 
there exist two different Wigner ctystals, with occupancy 
schemes ...1100... and ...1010..M At first glance, the 
extreme configuration ...1010..., corresponding to a pe- 
riod 2 "4kj?"CDWE3, appears to be strongly preferred, 
but in fact more careful analysis shows that it dominates 
the ground state only for fairly substantial VEil. This 
can be seen rigorously for U — ► oo, where the 1/4-filled 
spinful band cap be mapped rigorously to the 1/2-fillcd 
spinless bandE^I, which in turn can be mapped (via a 
Jordan- Wigner transformation) to an anisotropic Heiscn- 
berg spin 1/2 chainEj. Using this approach, one finds 
that the period 2 "4kp" CDW becomes the ground state 
only for V > Ve-r= 2 (in units of |t|)E3. For finite U, nu- 
merical resultsEj show that V c is slightly larger than 2. 
Given the estimated values of V in the organic CTS, it 
seems unlikely that they will exhibit this (...1010...) in- 
trachain ordering. This expectation is strongly supported 
by the result that the ...1010... CDW cannot coexist with 
the BOWtrD, whereas the (TMTTF) 2 X are known to ex- 
hibit.^, low-temperature transition to a SP-BOW ground 
stateM 

For V < V c the extended.. ID Hubbard model at 1/4- 
filling is a Luttinger liquidE^I that is also susceptible to a 
2kp bond and charge distortion, and it is this distortion 
that can be described by aay one of the four equiva- 
lent configurations ...1100... ED. The 2kp CDW compati- 
ble with the ...1100... configuration has the form 

Pc(j) = 0.5 + po cos(2kpja — 3ir/4) 

= 0.5 + p cos(7rj/2-37r/4), (5) 

This particular CDW also coexists with a BOW, since 
the charge-transfer across a '1 - 1' bond is different from 
that across a '1 - 0' (or '0 - 1') bond, which again is dif- 
ferent from the charge-transfer across a '0 - 0' bond. It is 
a subtle-,but crucial fact, confirmed by earlier numerical 
studiescZI, that this same CDW can now promote two dif- 
ferent BOWs, each with three different bond strengths. 
In each of these the '0 - 0' bond is the weakest, but de- 
pending upon the strength of the Coulomb interaction, 
the '1-1' bond can be stronger than a '1 - 0' (or '0 
- 1') bond (since charge-transfer in the former can oc- 
cur in both directions), but it can also be weaker (since 
charge-transfer in the former leads to double occupancy, 
while no double occupancy is created in the charge trans- 
fer between a '1' and a '0'). Consistent with this and the 
numerical resultsc3, we shall refer to the first bonding 
pattern as "SUWU" (for a strong '1-1' bond, undis- 
torted '1 -0' bond, weak '0 - 0' bond, followed by an 
undistorted '0-1' bond), where a strong bond has t$ > t, 
an undistorted bond has t\j = t, and a weak bond has 
t w < t. This BOW has pure period 4 "2kF" periodicity 
and is accompanied by lattice distortion 



Uj = uq cos(2kp ja — it/4) = uo cos(7rj/2 — tt/4). (6) 

Again consistent with the numerical results, we call the 
second bonding pattern "W'SWS" (for a stronger weak 
'1-1' bond, strong '1-0' bond, weak '0 - 0' bond and 
strong '0-1' bond, with t$ > t > tw> > tw)- Interest- 
ingly, the W'SWS bonding pattern is a superposition of 
the pure 2k^ period 4 SUWU structure and the pure 4kp 
period 2 SWSW structure and is accompanied by lattice 
distortion 

Uj = u [r 2 k F cos(2k F ja - tt/4) + r 4kF cos(4k F ja)} (7) 
= u [r 2 k F cos(7r//2 - 7r/4) + r ikp cos(tt/)], 

where r 2 k F and rik F are the relative weights of the 2k F 
and 4ki? bond distortions, respectivelycJ. These results 
were established numerically in reference [27], where from 
comparisons to available experimental data in the 1:2 an- 
ionic TCNQ systems it was also shown that the phase 
relationship between the coexisting 2k F CDW and the 
W'SWS BOW (the W bond connects sites with greater 
charge densities than the W bond) is precisely in agree- 
ment with theory. 

Very importantly, we show below that the dimerization 
of the dimer lattice with one electron per dimer also leads 
to a W'SWS bonding pattern (see Fig. 0(c)), which in its 
turn promotes the site occupancy scheme ...1100.... This 
coexistence will therefore occur in either the full 1/4- 
filled band model or the effective 1/2-filled, dimerized 
dimer approach. 

C. 1/4-filled band, t± « t, U,V ^ 0. 

The above two BOW-CDWs describe the ground state 
of the interacting 1/4-filled band in the limit of t± = 
0, where the '1 - 1' bond is a singlet. As in the 1/2- 
filled band though, singlets are expected to give way to 
SDW order for t± ^ 0. Thus we must understand the 
role of the spin degrees of freedom. Once specific spins 
are assigned to the sites labeled '1' in the ...1100... con- 
figuration, the sites labeled '0' become distinguishable, 
as a given '0' site is now .closer to one particular '1' (up 
or down) than the otherc.1 In this case the '0' site is 
expected to acquire the spin characteristic of its neigh- 
boring T'. The charge and spin along a chain can now 
thus be denoted as tilji.T, where the sizes of the arrows 
are schematic measures of the charge and spin densities 
on the sites. Note that this represents the SDW of the 
form 

Ps(j) = {c] M a c j-MA ~ c j,M,l c i,M,l) ( 8 ) 
= Ps2k F cos(2k F ja - it/ 4) + p sikp cos(4k F ja - %), 

which coexists with the BOW and CDW. 

Commensurability effects imply that the possible 
phase shifts between adjacent chains in the anisotropic 
2D system are 0, n/2 and it, and we have performed 
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FIG. 2. Sketches of the BCSDW ground states that oc- 
cur for small t± in the strongly correlated, anisotropic 2D 
1/4-filled band. The arrows indicate the spin directions and 
their sizes indicate the relative charge and spin densities. The 
hopping integrals used to calculate the energies of the dis- 
torted lattices correspond to (a) r^u F = (see text, Section 
V) and (b) r^u F = T2k F , and are shown above the bonds 
along the top chain. This variation in t reflects the BOW. 
The bond-distortion pattern in (b), with slightly modified 
weak bond hopping integrals, also corresponds to the dimer- 
ized dimer lattice for small enough t±. Note that the charge 
ordering corresponds to the ID paired electron crystal along 
the longitudinal and both diagonal directions and the the 
monatomic Wigner crystal along the transverse direction. 

explicit numerical calculations to determine that the 
lowest energy state is obtained with a phase shift of n. 
The intrachain bond orders, determined by the probabil- 
ities of nearest neighbor charge-transfers, continue to be 
different for the different pairs of neighboring sites. This 
is the major difference between the possible broken sym- 
metries in the 1/2-filled and 1/4-filled band. While in the 
1/2-filled band there is no overlap between the extreme 
configurations favoring the BOW, CDW and SDW, in 
the weakly 2D 1 /4- filled band the same extreme configu- 
ration supports all three broken symmetries^. For small 
nonzero t± , we therefore expect a strong cooperative co- 
existence between the BOW, the CDW and the SDW. 
Furthermore, since the same CDW coexists with both the 
SUWU BOW and the W'SWS BOW, this coexistence is 
independent of which particular BOW dominates. This 
has been explicitly demonstrated in reference [24] , where 
it was shown that the overall ground state for small t± 
is one of the two BCSDW states shown in Fig. |[ with 
overall 2D periodicity of (2kp, tt). 



D. 1/4-filled band, t±<t,U,V=£ 0. 

What happens as t± is further increased? Within k- 
space single-particle theory, increasing t± should destroy 
the nesting of the Fermi surface. But as we have indicated 
above, our real space analysis predicts, and our numer- 
ical results will establish, that this destruction does not 
occur. To argue this convincingly, we must first show 
how this destruction of the nesting, which certainly does 
occur for non-interacting electrons, can be correctly de- 
scribed within our configuration space picture of the bro- 
ken symmetry. Recall that the one-electron hopping term 
in Eq. (1) introduces "paths" between the extreme con- 
figurations, where each step in a given JtaJth connects two 
configurations related by a single hopoBo. Nonzero tj_ 
introduces many additional paths connecting the extreme 
configurations that are the 2D equivalents of ...1100... 
(with a 7r-phase shift between consecutive chains). For 
U = V — 0, there is no inhibition of these paths, and 
it therefore becomes easier to reach one extreme config- 
uration from another, leading to enhanced configuration 
mixing (relative to ID), which in its turn destroys the 
"nesting" and the broken symmetry. 

The situation described above changes, however, for 
nonzero Coulomb interaction. Interchain hopping t± 
leads to partial double occupancy on a single site (t |) 
with an energy barrier that, while less than the bare U, 
is a U e ff that increases with U. The energy barrier to in- 
terchain hopping leads to "confinement" of the electrons 
to single chains, a concept that has been widely de^ 
recently, in the context of high T c superconductors 
For large enough U e ff, the confinement can be strong 
enough that the broken symmetry state can persist up 
to the isotropic limit t± ~ t. 

More precisely, the bond and charge components of the 
BCSDW can persist up to the isotropic limit t± ~ t, lead- 
ing to the BCDW state we have previously introduced. 
The evolution of the spin structure is different from and 
more subtle than the bond and charge components. From 
the cartoons in Fig. we see that for the SDW to exist 
it is essential that the '0's have a spin "direction" . In 
the small t± case, the sign of the spin on a '0' is neces- 
sarily that of the nearest inteachain '1'. Note, however, 
that each '0' also has two interchain 'l's as neighbors and 
that for a stable SDW, the spin densities of the 'l's that 
are neighbors of a specific '0' must be opposite (as shown 
in the Figure). Therefore, with increasing tx, competing 
effects occur. On the one hand, the magnitude of the in- 
terchain exchange coupling J± ~ t\/U e ff increases. On 
the other hand, the spin density on a site labeled '0' de- 
creases because of the canceling effects of the intra- and 
inter-chain neighboring 'l's. We thus expect the SDW of 
the 2D lattice to vanish at a t\ that will depend on the 
magnitudes of the bare U and V. 

This description of the evolution of the SDW applies 
to the true 1 / 4-filled band. In lattices that are dimerized 
initially, further dimerization leads to the occupancy TO' 
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FIG. 3. Numerical results of ID simulations for [7 = 6, 
V = 1. (a) Charge densities (numbers inside each circle, 
which represents one molecular site) and bond orders (num- 
bers against the bonds) at the center of an open uniform chain 
of 16 sites for a = j3 = 0. (b) Bond orders in a 16 site peri- 
odic ring with uniform hopping, and with externally imposed 
period 4 magnetic field of the same form as in Fig. ^| with 
amplitude e = 0.05. (c) Same as in (b) with e = 0.1. Be- 
cause of equal bond lengths and nonzero V , there is a weak 
contribution by the ...1010... CDW to the ground state here 
and the charge densities are not pure ...1100... The filled 
(unfilled) circles correspond to large (small) charge densities. 
The bond orders also show weak deviation from pure SUSU 
or W'SWS behavior, and the bond orders shown are averages 
for each kind of bond. The magnetic field induced SDW cre- 
ates a spontaneous BCDW. (d) Charge densities in a periodic 
dimerized dimer lattice of 16 sites. The double bond corre- 
sponds to t = 1.2, and the dotted and double dotted bonds 
to t — 0.7 and 0.9, respectively. Note that the CDW pat- 
tern in this effective 1/2-filled band system is the same as the 
1/4-filled band lattices in (a), (b) and (c). 

or '01' on each dimcr. If the original dimerization 
is very strong, the spin on a given '0' will continue to 
be strongly influenced by the spin on its partner in the 
dimer, and t c ± at which the SDW vanishes in this case 
will be larger. 

The robustness of the BCSDW and the BCDW rela- 
tive to the uniform metallic state can be understood from 
the cartoon occupancy schemes in Fig. ||. It is instruc- 
tive to discuss the BCDW state in terms of the two large 
U Wigner crystal structures discussed by HubbardLj. 
We refer to the ...1100... electron arrangement as that 
of a "paired electron crystal", and the ..1010... as the 
"monatomic Wigner crystal." For the 3D_Ipw den- 
sity electron gas, Moulopoulos and AshcroftEZl showed 
that there exists an intermediate density range where 
the paired electron crystal has lower energy than the 
monatomic Wigner crystal, and the region < V < V c in 
our discrete lattice case can be thought of as intermedi- 
ate between the V = and V > V c . A striking feature of 
the BCSDW and the BCDW occupancy scheme is that 
it is a paired electron crystal along the chains (...1100..., 
periodicity 2kp ), a monatomic Wigner crystal transverse 
to the chains (...1010...., periodicity 4kp), as well as a 



paired electron crystal along both diagonals (...1100..., 
periodicity 2kp ). It is thus possible to predict that even 
in the presence of interactions not explicitly included in 
Eq. (1), the BCDW continues to persist. For instance, 
by enhancing the 4k^ charge ordering along the trans- 
verse direction, the nearest neighbor interchain Coulomb 
interaction V± will further enhance the stability of the 
BCDW. Similarly, the diagonal ...1100... charge order- 
ing implies that even the additions of hopping tdiag and 
Coulomb repulsion Vdiag along the diagonals will not de- 
stroy the BCDW state for realistic parameters: in par- 
ticular, Vdiag stabilizes the BCDW relative to the other 
Wigner crystal (...1010...) along both x and y directions. 

In the above our goal has been to predict a novel semi- 
conducting state that is more stable than the metallic 
state. Even if this semiconducting state is assumed, 
however, there is an additional surprise in our claim, 
viz., the dominance of the singlet BOW over the SDW 
for strong two-dimensionality in the interacting quarter- 
filled band. This is exactly opposite to what is observed in 
the 1/2-filled band. While in the half- filled band a single 
singlet-to-antiferromagnet transition occurs with increas- 
ing t± , for the 1 / 4-filled band, a second antiferromagnet- 
to-singlet transition is predicted at large t±. Since a full 
discussion of this second transition at this junction would 
interrupt the flow of the narrative, we defer it to Ap- 
pendix 1, which presents arguments based on variational 
concepts and valence bond theory to motivate this result. 

IV. NUMERICAL RESULTS 
A. Results for ID lattices 

Computational limitations will compel us to use fairly 
small lattices in 2D and will prevent us from studying 
dynamical phonons (even at a classical, self-consistent 
level). As a consequence, we will have to work with ex- 
plicitly distorted lattices, rather than allowing the distor- 
tions to arise naturally, as they would in larger lattices 
calculated with dynamical phonons. To provide justifica- 
tion for this approach, in this section we (a) extend-pur 
previous ID results obtained with nonzero a and /5cJ to 
zero e-ph couplings, to demonstrate that these bond and 
charge distortions are unconditional, and (b) show that 
the dimerization of the dimer lattice (see Eq. (2)) leads 
to the same CDW as the monatomic 1/4-filled band. 

It is known that in a sufficiently long open chain the 
bond orders and the charge densities at the center of the 
chain show the behavior in the long chain limit, even in 
the absence of the e-ph coupling. In Fig. ||(a) we show the 
exact nearest neighbor bond orders and charge densities 
at the center of an open undistorted chain of 16 atoms 
with all hopping integrals equal, for U = 6, V = 1. Note 
that both the BOW and the CDW show the 2k F mod- 
ulations discussed in section III, and appear in spite of 
uniform hopping integrals. 
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Second, we recall that in a purely ID system, a LRO 
SDW can occur only if an external staggered magnetic 
field is applied. We therefore incorporate an additional 
(external field-like) term 

Hsdw = - £[nj,lcos(2kpj) + n^icos(2k F j + tt/2)} 
j 

(9) 

and consider H + Hsdw for the 1/4- filled band with am- 
plitude e = 0.1. In reference [27] the same Hamiltonian 
was investigated for the case of finite bond distortion. 
Figs. ||(b) and (c) show the bond orders and CDW for a 
periodic ring (zero e-ph coupling and undistorted hopping 
integrals) with the SDW tin superimposed on it. Note 
that because of the periodicity, the bond orders are uni- 
form for the finite ring for e = 0. For e = 0.05 (Fig. |3|(b)) 
and 0.1 (Fig. ||(c)), the externally imposed SDW creates 
spontaneous BOWs with i4k F = and i4k F ^ 0, respec- 
tively. 

In Fig. [|(d) we show the charge densities on a peri- 
odic ring of 16 sites, now for the dimerized dimer lattice 
(the hopping integrals here are 1.2, 0.9, 1.2 and 0.7). 
The charge modulations (which appear entirely due to 
modulations of the interdimer bond orders) on the sites 
are exactly as in Figs. |](a)-(c), with the larger charges 
occurring on the sites connected by the stronger weak 
bond (the W bond, with tw> = 0.9). In discussions of 
the spin-Peierls transition within the effective 1/2-filled 
band (corresponding to the dimer lattice) , it is usually as- 
sumed that the electronic populations within each dimer 
cell remains uniform in the spin-Peierls state. Fig. ||(d) 
clearly shows that this is not true. 

B. Results for 2D lattices 

To confirm the expectations based on the qualitative 
arguments of Section III, we use exact diagonalizatioH. 
and Constrained Path quantum Monte Carlo (CPMC)E3 
numerical techniques to calculate for representative finite 
2D lattices: (i) the electronic energy gained upon bond 
distortion, 

AE = £7(0) - E(u jtM ), (10) 

where E(v,j m) is the electronic energy per site with fixed 
distortion w^m along the chains; (ii) the site charge densi- 
ties pj t M for the bond-distorted lattices; due to the coex- 
istence of the BOW and the CDW, measuring the CDW 
amplitude that results as a consequence of the external 
modulation of the hopping integrals is exactly equivalent 
to the measurement of the bond order differences in the 
charge-modulated lattices; and (iii) the z-z component 
of the spin-spin correlations, for a range of U , V and 
t_i_. We consider three distinct distorted lattices, two of 
which correspond to those shown in Figs. ||(a) and (b), 
where we have indicated the hopping integrals along the 



chain (the uniform lattice has a hopping integral of 1.0 
corresponding to all intrachain bonds). The third dis- 
torted lattice we consider is the dimerized dimer lattice, 
the hopping integrals for which will be discussed later. 

Ideally, calculations that aim to demonstrate persis- 
tence of a spatial broken symmetry should do fully self- 
consistent calculations of the total energy, which is a sum 
of the the electronic energy gain A£7 (including effects 
of both e-e and e-ph interactions) and the loss in lat- 
tice distortion energy. Unfortunately, in true many-body 
simulations (such as exact diagonalizations or CPMC) of 
the very large 2D lattices we investigate (see below) , such 
self-consistent calculations are not possible. A well-tested 
alternate approach! is to calculate only the electronic 
energy gain for fixed lattice distortion and compare the 
calculated AE against a known reference configuration, 
where the distortion is known to occur. This approach 
works because for a fixed distortion, the contribution of 
the elastic energy to the total energy is constant, inde- 
pendent of the other parameters; therefore the gain in 
electronic energy, relative to that for the reference config- 
uration, is a direct measure of the tendency to distortion. 
An example of a previous successful application of this 
approach is the enhancement by e-e interactions of the 
bond alternation in the ID 1/2-filled band; here, the ref- 
erence configuration corresponds to the limit of zero e-e 
interaction (SSH model), where the Peierls bond alterna- 
tion is known to occurH. For nonzero e-e interaction, the 
electronic energy gain for fixed bond alternation can be 
larger (see Figs. 2.26 and 2.31 in reference [j]) , indicating 
the enhancement of the bond alternation by e-e interac- 
tion, a theoretical result that has been confirmed by all 
subsequent studies. Similarly, in the 2D 1/2-filled band, 
calculations of the electronic energy gain for fixed bond 
distortion have been used to prove the decrease in the 
tendency to Peierls bond alternation upon the inclusion 
of e-e interaction (see Fig. 10 in referetige [To]), a result 
that is in agreement with other studieaS'EI as-well as the 
determination of long range AFM in this cased. Thus the 
approach has been shown to work in two cases in which 
exactly opposite outcomes, - in one case, an increase in 
dimerization, in the other decrease, occurred, in- 

dicating its robustness. 

At first glance, it appears that there exist two differ- 
ent reference configurations in the present case. First, for 
given t±, one could study A£7 as a function of U and V: 
in essence, this amounts to comparing uncorrelated and 
correlated lattices for each t±. Second, for given U and 
V, one could calculate AE as a function of t±. In fact, 
the first approach does not yield correct results for two 
reasons: (i) the uncorrelated 2D lattices are undistorted, 
so there is no obvious AE with which to compare the 
correlated results; and (ii) magnitude of AE decreases 
with U and V even in the ID limit, where we know that 
the bond and c harg e distortions are unconditional (see 
references |l6|- |l9| , p7| as well as the immediately previous 
subsection on ID numerical results). Thus to determine 
properly the tendency to distortion in 2D, our reference 
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configuration should be the single chain. We therefore 
normalize the energy gained for coupled chains (A.E) 
against that for the single chain (AEq) with the same 
U and V. A decreasing AE/ AEq as a function of t± sig- 
nals the destruction of the distortion by increasing two- 
dimensionality, while a constant, ar increasing AE/AEq 
indicates a persistent distortiontl'El. Since the BOW and 
the CDW are coupled cooperatively, the behavior of the 
charge ordering gives a second measure for the tendency 
to bond distortion. Decreasing charge ordering for fixed 
bond distortion, as a function of t± (as occurs for nonin- 
teracting electrons), indicates the tendency to decreasing 
bond distortion, while constant or increasing charge or- 
dering indicates persistent bond distortion. The expected 
(and calculated, see below) charge ordering pattern is the 
same for all bond distortion patterns and is the same as in 
ID (with, however, a 7r-phase shift between consecutive 
chains). 

As mentioned above, our numerical calculations in- 
volve both exact diagonalization and the CPMC tech- 
nique. Because of the sign errors that plague quantum 
Monte Carlo calculations in 2D, it is critical to obtain 
a precise idea about the accuracy of the numerical re- 
sults. This is especially so bepapae CPMC calculations 
that have been reported so farE3Ca are only for the sim- 
ple Hubbard Hamiltonian and did not include the near- 
est neighbor interaction V. In Appendix 2 we discuss 
our methodology and give detailed comparisons of ener- 
gies and correlation functions obtained for finite lattices 
within the CPMC and exact diagonalization procedures. 
As shown there, although the CPMC technique is not 
variational, the accuracies in both energy and correla- 
tion functions are sufficient for our purposes. 



8x2 U=6 V.1 

□ 8x2 U.6 V.1 
• 8x6 U.6 V.1 

O 16x6 U.6 V.1 

8X2U.V.0 

8x6 U.V.0 

16X6U.V.0 


I I 

I 




i — — <; 

i 

* 





1 1 I I I I I I I I 

0.2 0.4 0.6 0.8 



FIG. 4. AE/AEo versus t± for a 2k F bond distortion 
(i~4k F = 0) for the 8 x 2, 8 x 6, and 16 x 6 lattices for 
U = V = and for U = 6, V = 1. For the 8 x 2 lattice both 
exact (solid line) and CPMC results are shown. Intrachain 
hopping integrals for the distorted lattices are as indicated in 
Fig. |(a). 



For numerical results obtained from finite-size calcu- 
lations to be relevant in the thermodynamic limit, it is 
essential to choose proper boundary conditions. In the 
present case, we choose lattices and boundary conditions 
based on the physical requirement that for noninteract- 
ing electrons any nonzero t± must destabilize the BCDW 
on that particular finite lattice. Details of the analysis 
that guided our choice of 2D lattices are also presented 
in Appendix 2. There we show N x M lattices 
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FIG. 5. (a): Site charge densities on one of the 6 chains 
in a 2kp bond-distorted 8x6 lattice, for t± = 0.2, U — 6, 
and V — 1. The line is meant as a guide to the eye. Note 
the expected "..1100.." structure discussed in the text, (b) 
Amplitude of the 2k F CDW for the 2k F bond-distorted 8x2 
(exact), 8x6, and 16x6 lattices. The ground state of the 
16x6 lattice is in the S = 1 subspace for t± > 0.6, and the 
CDW amplitudes for the S = states here are expected to be 
greater than those calculated for the ground state and shown 
in the Fig. (see text). 

(with N the number of sites per chain and M the num- 
ber of chains) that obey the above physical requirement 
are restricted to those for which N = 8n, where n is an 
integer. On the other hand, there is no restriction on 
M, except that M be even to avoid even/odd effects. In 
our calculations below, we have chosen M = 4n + 2, for 
reasons that are also discussed in Appendix 2. 

We make one final point before presenting the 2D nu- 
merical data. The restriction to N = 8n sites coupled 
with the 1/4-filling introduces a potential subtlety into 
the numerical computations of AE/AEq for nonzero U 
and V. Finite 4n-electron non-l/2-filled ID undistorted 
periodic rings have their ground state in the total spin 
S = 1 subspace, and even the distorted system's ground 
state can be in the S = 1 subspace for the smallest 4n- 
electron rings. We have confirmed from exact 
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FIG. 6. The z-z spin correlations between sites 1 (left panels) and 2 (right panels) on the first chain of the 8x6 lattice and 
sites j = 1 - 8 on the second chain, with U = 6, V = 1 for four values of t±. Due to finite size effects the wavefunction has 
small admixing with the ...1010... charge order which affects the individual magnitudes of the spin-spin correlations (see text). 
AFM correlations increase with t± up to t± = 0.4 but then vanish at t± ~ 0.6, even though the BCDW continues to persist 
for all t± (see Fig. fj). Lines are guides to the eye. 



diagonalizations of the 8x2 lattice that the ground 
state is in the S = state for the smallest nonzero t±. 
Thus while AEq can correspond to the energy gained 
upon distortion in the S = 1 subspace, AE necessarily 
corresponds to the energy gained upon distortion in the S 
= subspace. As this important but subtle point requires 
extensive discussion that would interrupt the presenta- 
tion here, we present the details in Appendix 3, where we 
show that despite this subtlety, the behavior of AE/ AEq 
nevertheless is a proper measure of the stability of the 
distorted state for nonzero 1 1 . 



1. Exact diagonalization and CPMC calculations, r4k F = 

In Fig. || we show the behavior of AE/AEq for the 
non-interacting and interacting (U = 6, V = 1) cases for 
three different lattices satisfying our boundary condition 
constraints. In all cases we measure the electronic energy 
gained upon 2kp SUWU bond distortion (corresponding 
to nearest neighbor hopping integrals t$ — 1.14, tjj — 
1.0, and tw = 0.86), relative to that of the undistorted 
state with equal hopping integrals. For the 8x2 lattice 
the calculations involved both exact diagonalization and 
the CPMC technique. The 8x2 results, taken together, 
then provide an estimate of the precision of the CPMC 
calculation. The exact diagonalization studies also con- 



firm that the system is in the total spin state S = for 
t± as small as 0.01 (see Appendix 3). 

The large scatter in the normalized AE at very large 
and very small tj_ may be due to the degeneracies in the 
non-interacting system at t± — ► and t± — ► 1. Further- 
more, as pointed out in Appendix 2 (subsection A), the 
absolute values of AE are rather small, especially for the 
pure Zkp (f4fc F = 0) distortion. The systematic errors 
due to the CPMC approximation are therefore large in 
these two regions. Nevertheless, except for the AE/AEq 
value at t± = 0.1 for the 8x6 lattice, at all other t± 
the AE/AEq values are above 1 for all three lattices, 
and far above the normalized non- interacting values. As 
seen in Fig. j|, while for the non-interacting cases the 
AE/AEq decreases rapidly with t±, for the interacting 
cases the AE/AEq either remains unchanged or is en- 
hanced by t ± . Because of the strong degeneracies in the 
one-electron occupancy scheme at the Fermi level at t±_ = 
1, a single well-defined one-electron wavefunction is miss- 
ing here. The CPMC calculations therefore could not be 
done for t± — 1.0. It is, however, highly unlikely that the 
BCDW persists for t± = 0.9 but vanishes at t± = 1; this 
expectation is corroborated by the results of the exact 
diagonalization studies for the 8x2 lattice, which were 
performed for the full range of t±, including tj_ = 1 and 
showed enhanced distortion throughout the whole region. 
In the following sections we also show AE/AEq for the 
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FIG. 7. The z-z spin correlations between sites 2 (left panels) and 3 (right panels) on the first chain of the 16 x 6 lattice and 

sites j = 1 -16 on the second chain, with U = 6, V = 1 for four different values of t±. The finite size effects, and contamination 

with the ...1010... charge order here is smaller than in Fig. ^| Lines are guides to the eye. 

2kp+4kp (r4fc F ^ 0) and dimerized dimer lattice. In the 16x6 lattice is a spin effect: the bond distorted state 
both of these cases, the magnitude of AE is larger and 
hence easier to compute, but degeneracies restrict CPMC 
simulations to smaller t±. In both cases, AE/AEo is 
close to or above 1 for all t± we have studied. 

As discussed in the above, the bond-distorted lattices 
(both r^kp = and r^k F ^ 0) have a synergetic coexis- 
tence with the CDW. Thus the amplitude of the CDW, 
defined as Ap c = p c i—p C s, where p c i and p cs are the larger 
and smaller charge densities on the ...1100... 2kp CDW, 
is an alternate measure of the stability of the BOW. If the 
nonzero t± destabilized the bond-distortion, then even 
with fixed 2kp distorted hopping integrals the amplitude 
of the BOW (measured as the differences in the bond 
orders) would decrease, and the diminished strength of 
the BOW in turn would decrease Ap c . This is easily 
confirmed for the noninteracting Hamiltonian, where the 
amplitude of the CDW decreases with increasing t ± . In 
Fig. |](a) we show the charge densities on a single chain 
for a bond-distorted 8x6 lattice (because of periodicity, 
all chains are equivalent) for U = 6, V = 1, and t± =0.2. 
In Fig. ||(b) we have shown the behavior of Ap c for all 
the three lattices we have studied, now as a function of 
t j_ . Degeneracies in the one-electron energy levels in the 
16x6 lattice for t± > 0.6 even with finite bond-distortion 
cause the CPMC ground states in this region to be S = 
1. Exact calculations in the ID limit show that the am- 
plitude of the CDW in S — 1 is less than that in S = 
0. Thus the weak decrease in the Ap c values with tj_ in 



is S = at small t± and S = 1 at large t±. The Ap c 
values at large t± for the 16 x 6 lattice should there- 
fore be considered as lower limits (the Ap c values of the 
16 x 6 lattice are considerably larger than that of the 
S = 1 single chain of 16 sites). In agreement with the 
behavior of the AE in the interacting case (see Fig. ||), 
the CDW amplitude now increases or remains constant 
with increasing t± for all the lattices studied, indicating 
a greater tendency to bond and charge distortion with 
increasing t±. Taken together, the results of Figs. and 
provide quantitative proof of our qualitative arguments 
establishing that the B 'CDW is a robust broken symmetry 
state for the interacting 2D j-filled band. 

In Fig. ^ we show the inter-chain spin-spin correla- 
tions between sites 1 and 2 on the first chain, and sites j 
= 1 - 8 on the second chain, for the 2kp bond-distorted 
8x6 lattice for several values of t±. The SDW profile 
is somewhat different from what is expected from a pure 
...1100... charge modulation along the chains because the 
wavefunction of this finite lattice also has contributions 
from the ...1010... type intrachain charge modulation. 
The small ...1010... contribution to the wavefunction af- 
fects the charge density, %,M,t + n j.M.i only weakly, but 
the spin density, being the difference nj,M,T ~ n j.M.i is a 
smaller quantity and is affected relatively more strongly. 
It is useful here to recall however that within the rectan- 
gular lattice, ...1010... charge orderings along both lon- 
gitudinal and transverse directions give triangular lattice 
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of occupied sites, and thus a pure ...lOlO-.p-cannot give 
the SDW profiles of Fig. | (see also below)E3. 

Qualitatively, at t± = 0.1 the SDW behavior is the 
same as in Ref. where these calculations were done 
for the 12x4 lattice: the amplitude of the interchain spin- 
spin correlation is independent of the distance between 
the sites, indicating long-range order. The qualitative 
behavior of the spin-spin correlations is the same for t± 
= 0.4, where, however, the amplitude of the SDW is 
larger. At still larger t±(= 0.6), the inter-chain corre- 
lations are very strongly antiferromagnetic at short dis- 
tances (j = 1,2 on chain 2), but the antiferromagnetic 
correlations have disappeared at larger distances. This 
can be seen from comparisons of the spin-spin correla- 
tions corresponding to values of j lying near the center 
of the second chain (j = 5), which are farthest from the 
spins occupying sites 1 and 2 on the first chain. While 
the spin-spin correlations near j=5 increase from t± = 
0.1 to 0.4, they decrease as t± is further increased to 0.6. 
Similarly, focusing on site 8 of the second chain, we see 
that the spin-spin correlation with site 1 on the first chain 
has actually changed sign upon increasing f j_ to 0.6 from 
0.4 (due to the very strong short-range antiferromagnetic 
correlations), and the magnitude of the positive spin-spin 
correlation with site 2 on the first chain has decreased. 
All of these results indicate the absence of long-range 
spin order for large t± > 0.6 in the 8x6 lattice. The 
loss of the long-range spin-order is most clear at t± = 
0.9, where spin-spin correlations are nonzero only for the 
nearest interchain neighbors. 

Fig. [7] shows the inter-chain spin-spin correlations be- 
tween sites 2 and 3 on the first chain and sites j=l. . .16 
on the second chain for the 16x6 lattice. The admix- 
ture of the intrachain ...1010... CDW is weaker in this 
larger system: this is because the "tunneling" between 
the extreme configurations ...1100... and, say, ...0110..., 
decreases with size, and as consequence, V c increases with 
size in finite systems. This can be seen by simply com- 
paring the figures on the left and right panels for t± — 
0.1 and 0.2. If the intrachain CDW were a pure ...1010..., 
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FIG. 8. AE/AEo versus t± for r ikp = r 2 k F for the 8 x 
2 and 8x6 lattices for U = 6, V = 1. Intra-chain hopping 
integrals for the distorted lattices are as indicated in Fig. 0(b) . 



the signs of the spin-spin correlations for each j would 
be the same for both i = 2 and i = 3. Different signs for 
these correlations are signatures of the ...1100... CDW 
(see Fig. ||). As in the 8x6 system, long-range SDW be- 
havior is seen for t± = 0.1. Focusing on sites j=7 - 12 
on the second chain, the amplitude of the SDW increases 
from t± = 0.1 to fj_ = 0.2, but further increasing t± 
to 0.3 destroys the long-range order, as evidenced again 
by very large AFM correlations at short distances and 
vanishing correlations at large distances (sites j=7. . .12 
on the second chain) . The vanishing of the SDW is seen 
most clearly at very large t± (tj_ = 0.9 in Fig. |?]). We 
observe this same behavior of the SDW on 8 x 2 lat- 
tice. In all cases, the SDW amplitude initially increases, 
exhibits a maximum, and then vanishes at a ij_ which 
decreases with the size of the system. As discussed in 
section III.D, this behavior is to be expected from the 
nature of the BCSDW in Fig. |. The initial increase of 
the SDW amplitude indicates that t\ is nonzero, a con- 
clusion that is also in agreement with the experimental 
observation of the BCSDW state in the weakly 2D or- 
ganic CTS (see below). Based on the calculations for 16 
x 6 lattice, we estimate 0.1 < t\ < 0.3 for the strictly 
rectangular lattice for U = 6, V = 1. 



2. Persistent distortions with r4k F / 

The bond modulation pattern in the 1/4-filled band 
given in Eq. (7) has in general both 2kp and Akp com- 
ponents. Figs. ^ and || show persistent distortion at 
large inter-chain couplings for = (purely 2kp bond 
distortion). The persistent BCDW is expected also for 
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FIG. 9. AE/AEo versus t± for a dimerized dimer lattice 
for the 8x2, 8x6 and 16x6 lattices, for U = 6, V = 1. The 
intradimer hopping integrals are 1.2 in both cases. All inter- 
dimer hopping integrals are 0.8 in the dimer lattice, and 0.7 
and 0.9 in the dimerized dimer lattice. 
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FIG. 10. The z-z spin correlations for the 16x6 dimerized dimer lattice. Correlations are shown between site 2 and site 3 on 
the first chain and sites 1-16 on the second chain with U — 6, V = 1 for two values of t± (0.1 and 0.37). Lines are guides to 
the eye. 



Physically, the reason for this persistence is the coesisi- 
ing site CDW, whose nature is independent of r^Ej'EH. 
We show in Fig. ||the calculated AE/ AE for r 4 fe F = r 2 k F 
(equal admixtures of 2kj? and 4kp bond distortions), for 
the 8 x 2 and 8 x 6 lattices for U = 6 and V = 1. The 
hopping integrals corresponding to the distorted lattice 
here are 1.089, 0.974, 1.089 and 0.848, and the energy 
gained is being measured against the uniform lattice. 
Starting from t± = 0.5, the one-electron AE is highly 
discontinuous. This is because distortions with r^p ^ 
do not correspond to a natural periodicity for the non- 
interacting system. As a consequence the noninteracting 
wavefunctions are not suitable trial wavefunctions for the 
CPMC calculation. For the same reason the 16x6 cal- 
culation could not be performed here. The similarities 
between the results for the 8x2 and the 8x6 lattices 
are obvious. The ratio AE/AEo is independent of t± over 
a broad range of t± and increases slightly for large t±, 
indicating once again a stable 2D BCDW. Although only 
limited data could be obtained for this case, the dimer- 
ized dimer lattice is very similar in character to /'4^ F =/= 
(see Fig. ||(b)). In the following we show convincing 
evidence for persistent double-dimerization in 2D. 



3. The dimerized dimer lattice 

We have previously noted that Fig. ||(b) suggests that 
an alternate way to view the BCDW/BCSDW states is 
as a dimer lattice with additional structure within each of 
the dimer cells; the dotted box in Fig. ^(b) represents one 
dimer. Each dimer has ona-electjcan, leading to an "effec- 
tive half- filled" dimer bandO ~c2e3. Bond dimerization. in 
the ID 1/2-filled band is unconditional for all U > 2VM, 
and thus this dimer lattice itself distorts in a period 2 
dimerization pattern in ID. In this section we show the 



additional result that the (anisotropic) 2D dimer lattice 
is unconditionally unstable to a second dimerization for 
all t±. 

We choose the hopping integrals between the two sites 
within the dimer cell to be 1.2 in our calculations. The 
two inter-dimer hopping integrals for the uniform dimer 
lattice were taken to be 0.8, while for the distorted 
("dimerized") dimer lattice these were taken to be 0.7 
and 0.9, respectively (i.e., the dimerized dimer lattice has 
hopping integrals 1.2, 0.7, 1.2, 0.9 along each chain). Ex- 
act diagonalizations show that a 7r-phase shift between 
the chains (i.e., dimer cells lying directly above each 
other, but a strong inter-dimer bond on one chain fac- 
ing a weak inter-dimer bond on the next chain) gives the 
lowest total energy. Again we define AE and AEo as 
the electronic energies gained per site upon intcrdimer 
bond distortion by the 2D and ID lattices. Fig. [9] shows 
the AE/AEo behavior for the 8x2 lattice over the com- 
plete range of t± and for the 8x6 and 16x6 lattices for 
several different t± for U — 6 and V = 1. The 8x6 
and 16x6 lattices, taken together, cover nearly the full 
range of t±, and the AE/AEo behavior for these lat- 
tices closely follow the curve for the 8x2 lattice. As 
before, AE/AEo is significantly greater than 1 for the 
complete range < t± < 1, indicating the persistence of 
the dimerization of the dimer lattice in the interacting 
case, whereas for the non-interacting case, the dimeriza- 
tion vanishes, as expected. 

Fig. [l^ shows the interchain spin-spin correlations be- 
tween sites 2 and 3 on one chain and sites j = 1 - 16 
on a neighboring chain, for a 16x6 dimerized dimer sys- 
tem. Notice the far smaller contribution by the ...1010... 
intrachain charge ordering here. This is because of the 
large difference between the hopping integrals even in 
the "uniform" lattice with interdimer hopping integrals 
of 0.8 here. Such a large bond dimerization diminishes 
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the intrachain ...1010... contribution. The spin-spin cor- 
relation amplitudes cannot be directly compared to Fig. |7 
because of the different distortion amplitudes, but Fig. [lC 
shows that the SDW amplitude is significantly greater in 
the intermediate t± regime (t± = 0.37 in the Figure) com- 
pared to the small t± regime unlike the results in Fig. |t]. 
Our calculations indicate that the larger the difference 
between the intra-dimer and the inter-dimer hopping in- 
tegrals, the greater the range of the t± over which the 
SDW is stable. Thus with hopping integrals of 1.2, 0.9, 
1.2 and 0.7 along each chain, the SDW in the 8x6 lattice 
persists even at t± = 0.6 (in contrast to the 2kp bond- 
distorted lattice of Fig. |2|) , but vanishes at still larger t±_ . 
This is expected from our discussion of the behavior of 
the SDW in Section III.D. Recall that the smaller spin 
densities on the sites labeled '0' are influenced by both 
the intrachain nearest neighbor as well as the interchain 
nearest neighbor with opposite spin, and this competi- 
tion creates a disordering effect. The larger the hopping 
integral between the '0' and the nearest intrachain '1', 
the larger the t± necessary to create the disordering of 
the spin, hence the greater stability of the SDW. We 
shall later argue that this same phenomenon is related to 
the very large magnetic moments of the k-(BEDT-TTF) 
salts. 



4- Effects of additional Coulomb interactions 

Fig. U clearly suggests that interchain nearest neighbor 
Coulomb interaction V± stabilizes the BCDW further. 
We have confirmed this by exact numerical calculations 
for the 8x2 lattice, as shown in Fig. |ll| below, where 
we have plotted AE/AEo for three different values of 
V±: 0, 0.5 and 1. Nonzero V± increases AE further. 
Similar calculations were done also with variable V± but 
fixed An even larger increase in AE is found in 

this case. Implementing V±_ over and above V is difficult 
within the CPMC, and therefore these calculations could 
not be performed for larger lattices. However, based on 
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FIG. 11. AE/AEo vs. t± for the 8 x 2 lattice, with U = 6, 
V = 1 and V± = 0, 0.5, 1.0, for r ikp = 0. 

the similarities between the AE behavior of the three 



lattices studied in Figs. || and [| no difference in the larger 
lattices is expected. 



V. SUMMARY OF THEORETICAL RESULTS 

We have performed detailed numerical calculations of 
various broken symmetries for the 2D 1/4- filled band 
within Eq. (1) and for the effective 1/2- filled band of 
dimer lattice within Eq. (2), for U = 6, V = 1. Re- 
garding these parameter values, the broken symmetries 
we have found will occur for all intermediate tcwstrong U 
but require V to be less than a critical V c > 2tEil. 

We have discovered three distinct new results in ID. 
First, we have confirmed that the BCDW state oc- 
curs spontaneously even for zero e-ph couplings (see 
Fig. ||(a)). The bond distortion pattern in the center of 
a long open chain corresponds to a pure 2k^? distortion, 
and coexists with the 2kp ...1100... type charge mod- 
ulation. Second, we have shown that a BOW appears 
spontaneously in a uniform periodic ring when the SDW 
t, is superimposed, confirming the synergetic coop- 
eration between e-e and e-ph interactions. The BOW 
pattern corresponds to = (see Eq. (7)) when the 
amplitude of the superimposed SDW is relatively weak 
(Fig. ||(b)), but switches over to r4k F ^ when the SDW 
amplitude is large (Fig. ||(c)). Our earlier demonstrations 
of the BOW-SDW coexistence were only for the bond dis- 
torted periodic systems. Finally, from exact calculations 
for a periodic dimcrized dimer ring, we have established 
the new result that the BOW here also coexists with the 
...1100... 2kp CDW, with the large (small) charges oc- 
cupying the sites connected by the stronger (weaker) in- 
terdimer W (W) bond (see Fig. ||(d)). Our earlier work 
had claimed that a 1/ 4- filled description was essential to 
obtain the BCDW and the BCSDW states. As shown 
in Fig. |3|(d), the same result is obtained, however, even 
for the dimer lattice, provided the second dimerization is 
allowed to occur. 

Three different bond distortion patterns were investi- 
gated in 2D. These correspond to r 4 fc F = (Fig. ||(a)), 
~£4k F = ?2k F (Fig. ||(b)), and the dimerized dimer lattice. 
In all cases a 7r-phase shift in the bond distortion be- 
tween consecutive chains gives the lowest energy. From 
calculations of energy gained upon bond distortion, we 
conclude that 2D bond distorted lattices with t4k F = 
and T4k F = ?2k F are both more stable than the uniform 
lattice (see numerical results in Figs. [| and ||). Simi- 
larly, the dimerization of the dimer lattice is also uncon- 
ditional (see numerical results in Fig. ||). The persis- 
tence of the distortions is a novel effect of e-e interac- 
tions and is in contradiction to what is expected within 
one-electron nesting concepts. The ground state of the 
strongly correlated 1/4-filled band is therefore a novel 
insulating BCDW state for all t±. 

The persistence of the BCDW for all anisotropies is 
also evident from the charge density calculations. In 
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Fig. ||, we have shown the amplitude of the CDW that 
accompanies the T4k F = BOW as a function of t±. In 
the absence of e-e interaction, the CDW amplitude de- 
creases rapidly with t± even with nonuniform hopping 
integrals. One interesting aspect of these calculations is 
that the CDW pattern is the same for all bond distor- 
tion patterns. Our computer capabilities do not allow us 
to determine self-consistently which of the three BOW 
patterns dominate within Eqs.(l) and (2) for a given U, 

V, tj_, a and (3. This is, however, largely irrelevant, be- 
cause the charge ordering is the same with all the bond 
distortion patterns. 

The SDW behavior is different from those of the BOW 
and the CDW. As seen from our numerical calculations 
of interchain spin-spin correlations in Figs. | and 0, the 
SDW amplitude of the novel BCSDW state is initially 
enhanced by t ± , but with further increase in t± the SDW 
vanishes, indicating a singlet BCDW state again in the 
large t± region. The range of t± within which a stable 
SDW is found depends on the BOW pattern, and within 
the dimerized dimer lattice (see Fig. M) the SDW can be 
stable over a wider range of t±. 

VI. COMPARISON TO EXPERIMENTS ON THE 
INSULATING STATES IN 2:1 ORGANIC CTS 

Experimentally, the organic cationic CTS, with 
catiomanion ratio of 2:1, span the range t± < 0.1 in 
(TMTTF) 2 X to t± ~ 1 in certain (BEDT-TTF) 2 X. 
Hence these materials provide a critical testing ground 
for our theoretical results. In reference [24], we com- 
pared our theoretical predictions regarding the BCSDW 
state to the mixed CDW-SDW found experimen- 
tally in (TMTTF) 2 Br, (TMTSF) 2 PF 6 and a-(BEDT- 
TTF) 2 KHg(SCN) 4 . Here we make additional, more de- 
tailed comparisons, distinguishing between ID TMTTF 
and weakly 2D TMTSF-based compounds, and also em- 
phasizing the similarities and differences between the 
salts of BEDT-TTF and BETS with different crystal 
structures. In the case of the TMTTF and TMTSF 
band structure calculations of hooping integrals have 
been been summarized by YamajEj. In both cases the 
lattice is anisotropic triangular in nature, which would 
correspond to our rectangular lattice with one additional 
diagonal hop tdiag beyond the usual t j_ . Both t± and tdiag 
are small in the ID TMTTF, while they are comparable 
in TMTSF and about 0.1 1£| in magnitude. As discussed 
in section III.D, the paired electron crystal ordering even 
along the diagonal directions in the configurations shown 
in Figs. 2(a) and (b) indicate that the BCDW and the 
BCSDW states continue to be stable for nonzero tdiag 
and there is thus no loss of generality in considering a 
rectangular lattice. Several crystal structures occur in 
the BEDT-TTF systems, and more subtle and individ- 
ual analyses for the different cases are required. Our aim 
is to show that a variety of recent experiments indicate 



that the BCSDW and the BCDW are appropriate de- 
scriptions of the insulating states of this entire class of 
2:1 cationic CTS, and conversely, the very nature of the 
insulating ground state in certain cases provides direct 
verification for some of our more surprising theoretical 
results. We discuss below each class of material individ- 
ually. 

A. (TMTTF) 2 X 

The (TMTTF) 2 X compounds are nearly ID semicon- 
ducting materials with weak to moderate dimerization 
along the stacks at high temperature. Because of this 
dimerization, they have often beeji-described within the 
effective 1/2-filled band pictureE3l23. Further dimeriza- 
tion of the dimerization occurs below the SP transition 
temperature Tsp (~ 15 K\_Existing theories of the SP 
transition in these systemso do not discuss the simul- 
taneous appearance of the 2kp CDW and assume that 
the site populations continue to be uniform below Tsp- 
As depicted in Fig. 0(c), and as confirmed in Fig. ||(d), 
independent of whether these systems arc considered 
as 1/4- filled or effective 1/2-filled with a dimer lattice, 
the appearance of this 2kp CDW is unconditional and 
the site populations are therefore not uniform. In a 
recent NMR study of 13 C spin-labeled (TMTTF) 2 PF 6 
and (|XMTTF) 2 AsF 6 charge-ordered states have been 
foundE3. Although such a charge-ordering suggests agree- 
ment with the theory presented here, one problem is that 
the initial appearance of the charge-ordered phase (at ~ 
70 K in_(TMTTF) 2 PF 6 ) occurs considerably above T SP 
(15 K)E3. There are two possible reasons why the charge- 
ordering might appear at a temperature Tco > Tsp- 
First, this might be due to fluctuation effects associated 
with the ID nature of the crystals. As has been shown 
by SchulzE-3, fluctuation effects associated with the SP 
transition may be seen at temperatures as high as ATsp, 
in which case signatures of charge ordering would also 
become visible at these high temperatures. The obser- 
vation of diffuse X, -ray . scattering at 2k P in this material 
already at ~ 60 seems to support this possibility. 

A second possibility is that the charge-ordering is driven 
primarily by the Holstein e-ph coupling (3 in Hamiltonian 
(1), and the SSH coupling a is small, such that actual lat- 
tice displacement and spin singlet formation takes place 
at lower temperature. Independent of which mechanism 
dominates to give Tco > T$p, it is important to keep 
in mind that (a) no charge-ordering is expected at all 
within conventional theories of SP transition, and (b) 
as discussed extensively in section III, charge ordering 
of the type ...1010..., as has sometimes been suggested 
(see below and footnote ^lj), promotes equal intrachain 
bonds, and therefore the SP transition could not occur if 
the ...1010... charge-ordering had takea-place. Finally as 
has been pointed out by us previouslyED, charge-ordering 
of the type ...1100... also occurs in the SP phase of the 
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anionic 1:2 TCNQ solids. 

Although most (TMTTF) 2 X exhibit the SP transi- 
tion, the material (TMTTF) 2 Br exhibits a transition 
to a SDWBEI, like the (TMTSF) 2 X. Also like the 
(TMTSF)2X, this material can become superconduct- 
ing, although at a relatively high pressure of 26 kbar. 
Within_the structural classification scheme described by 
JeromeE], this difference is due to the larger t± in 
(TMTTF) 2 Br (relative to the other TMTTF). We there- 
fore discuss this material along with the (TMTSF) 2 X. 



B. (TMTTF) 2 Br and (TMTSF) 2 X 

X-ray scattering studies by Ravy and Pouge t@0 have 
shown that in both (TMTTF) 2 Br and the prototype 
TMTSF system, (TMTSF) 2 PF 6 , CDW distortions occur 
below the SDW transition temperature Tsdw- Similar 
conclusions have been reached also by Kagoshima et al.Cj. 
In (TMTTF) 2 Br evidence for a 4kp lattice instability was 
foundEElo, clearly suggesting that the insulating state 
here is the BCSDW of Fig. |(b). In (TMTSF) 2 PF 6 the 
authors claim a "purely electronic CDW", which would 
indicate the dominance of the 2kp CDW over the BOW. 
Since, however, in both the 1/4- filled band and the effec- 
tive 1/2-filled band, the 2k P CDW necessarily coexists 
with a BOW, the experimental work merely indicates 
that the transition to the BCSDW state is driven mainly 
by the Holstein e-ph coupling in Eq. (1) rather than the 
SSH coupling (i.e., a is small), so that the actuaLmodu- 
lations of the intermolecular distances are smallEJ. This 
would agree with one of the two possible reasonings given 
by us for T co being larger than T SP in (TMTTF) 2 PF 6 
and (TMTTF) 2 AsF 6 , as discussed above. 




FIG. 12. Schematic view of-ptructure of a-(BEDT-TTg) 
donor plane from Mori et al.LD and Ducasse and Fritschu. 
The solid lines correspond to stronger interstack t p hopping 
integrals, the dotted lines to weaker intrastack t c hopping 
integrals. The a and c directions indicated are the crystal 
axes, and the x and y directions correspond to along the chain 
and perpendicular to the chains in Fig. |2| The arrangement 
of the spins in the BCSDW state is indicated. Any SDW 
should be weak because of the nearly isotropic 2D nature of 
the lattice, but nonvanishing because of the nonzero t c , which 
becomes tdiag in the x-y coordinate system of Fig. 2 (see text). 



One additional comment appears to be necessary. 
Frohlich mode sliding conductivity has been seen in 
(TMTSF) 2 XE3. While this indicates a weak incommen- 
surability of the density wave (sec below), an equally im- 
portant point is that the sliding conductivity in the past 
has been ascribed to a SDW: the SDW collective trans- 
port is viewed as that of two CDWs, one for each spin 
subband. The actual displacement of the charge density 
is difficult to visualize in configuration space within this 
picture. We believe that the experimental demonstration 
of the coexisting CDW and the present theoretical work, 
taken together, suggest the more coherent viewpoint that 
the sliding mode conductivity is that of a BCSDW. 



C. Q-(BEDT-TTF) 2 MHg(SCN) 4 

This class of materials, with M = K, Rb, Tl and 
NH4 has been of considerable interest recently. M = 
NH4 is a superconductor, but M = K, Rb, Tl are non- 
superconducting. Early magnetic susceptibility studies 
in the M = K material had indicated anisotropic sus- 
ceptibility below the so-called "kink" transition that oc- 
curs at 10 K, indicating a SDW; here the kink refers 
to the change in slope that occurs in the tempera- 
ture dependence of the resistivity and the Hall coef- 
ficient. On the other hand, analysis of the angle- 
dependent magnetoresistance oscillations by Sasaki and 
Toyota led these authors to conclude already in 1995, 
prior to the experiments by Pouget and Ravy in the 
(TMTSF) 2 PF 6 , that the dominant broken symmetry in 
a-(BEDT-TTF) 2 MHg(SCN) 4 is a CDWl£i Since, how- 
ever, a CDW would not explain the anisotropic suscepti- 
bility, Sasaki and Toyota concluded that the broken sym- 
metry here is a "mysterious" state that is a "SDW accom- 
panied by a CDW" or a "CDW accompanied by a SDW" . 
Muon spin resonance studies indicate very small mag- 
netic, moment per BEDT-TTF molecule here, ~ 0.003, 
/i S C3 (to be compared against 0.08 pL B in (TMTSF) 2 XEl 
and 0.4 - 1 /i B per BEDT-TTF dimer in k-(BEDT- 
TTF) 2 Cu(CN) 2 ClEl, see below). More recent 13 C-NMR 
studies in the M = Rb indicate even smaller magnetic 



moment (if it .exists at all) ~ 1 10 4 /i_bL3. Re- 



cent theoreticatfl and experimental^ investigations con- 
clude cither that the dominant broken symmetry here is 
a CDW or that it is not a conventional SDWO. 

We point out here that a mixed state with very small 
magnetic moments is exactly what is expected within 
our theory. In Fig. [l2] we have given a schematic 
view of the structure of the donor plane in a-(BEDT- 
TTF) 2 MHg(SCN) 4 . The one-electron hopping integrals 
(called "t p " and "t c " in the figure) have been calcu- 
lated-using approximate one-electtan techniques by Mori 
et al.L3 and Ducasse and FritscbEa. Here the t p corre- 
spond to the interstack hopping and the t c to the in- 
trastack hopping. Four slightly different p-type integrals 
and three slightly different c-type integrals are obtained 
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by these authors. We ignore the small differences within 
each type of hopping integrals, as a more important effect 
is the periodic modulation that appears with the BCDW. 
We believe that what is relevant in the present context 
is that t p > t c . The a-BEDT-TTF lattice is then simply 
a rotated (by approximately 45°) version of our rectan- 
gular lattice with both t and t± = t p and tdiag = t c . Our 
calculations (see Figs. |[ [| and ^) show that even 
at t± ~ 1 the correlated 1 /4-filled band (or the dimer- 
ized dimer lattice) remains bond and charge-distorted, 
while based on the ..1100... ordering along the diagonals 
we have argued that tdiag does not destroy this order 
(see section III.D). Furthermore, while t± > t\ destroys 
the SDW order (leaving the BCDW intact) by disorder- 
ing the spins on the sites labeled '0' (see section III), 
a small tdiag will have a tendency to restore it, since 
now each small spin has two neighbors with spins of the 
same sign and one spin with opposite sign. Thus, the ex- 
perimentally observed strong BCDW and a weak nearly 
vanishing SDW is exactly what we expect within our the- 
ory. Further evidence for a partial gap has been found 
in the 13 C-NMR studies of a-(BEDT-TTF) 2 KHg(SCN) 4 
in high magnetic fields, in a regiop-where the system was 
previously thought to be a metal_. In Fig. |l2] we give a 
schematic of the spin arrangement in the a-BEDT-TTF 
lattice; note that the underlying x «-+ y symmetry in the 
isotropic 2D limit implies that there are two degenerate 
orthogonal 2D BCDW states here. 

Since in a-(BEDT-TTF) 2 MHg(SCN) 4 charge-oxdering 
has also been discussed, by Kino and Fukuyamao, and 
more recently, by Sea_, we should point out that the 
charge-ordering proposed by these authors is different 
from that in Fig. Our charge-ordering in Fig. 1^ is 
a rotated version of Fig. || where the occupancy scheme 
is ...1100... along the x-direction and along the diago- 
nals. The charge-ordering found by Kino and Fukuyama, 
and by Seo, assumes that the ...1010... order dominates 
over the ...1100... order. The ordering determined by 
Kino and Fukuyama is within a Hartree-Fock solution to 
the simple Hubbard model (zero intersite Coulomb in- 
teraction and zero e-ph coupling) and consists of a stripe 
structure with stack occupancies (c-direction in Fig. |l2] ) 
alternating, i.e., stacks are either completely filled or 
completely devoid of holes). More recently, Seo has re- 
peated these calculations by incorporating nearest neigh- 
bor Coulomb interaction V, but by treating U within 
the Hartree-Fock approximation and the V within the 
Hartree approximation. Different stripe structures, in- 
cluding that of Fukuyama and Kino, are found now, but 
once again, these are derived fundamentally from the oc- 
cupancy scheme ...1010.i_.As has, however, been pointed 
out by previous authora_'_l, the ...1010... charge order- 
ing for the case of V — is an artifact of the Hartree- 
Fock approximation. Similarly, the Hartree approxima- 
tion for V also exaggerates the ...1010... order while the 
Hartree-Fock treatment of the Hubbard term exagger- 
ates the SDW order_l This is precisely why these au- 
thors find very large magnetic moments in the a-phase 



materials, in disagreement with experiments. 



D. k-(BEDT-TTF) 2 X 

The deviation from the rectangular lattice is much 
stronger heret_. Crystal structure effects are very strong, 
and as a consequence the lattice is strongly dimer- 
ized, with-the dimer sites forming an effective triangu- 
lar latticet_. The strong deviation from the rectangu- 
lar lattice precludes direct comparisons against our the- 
ory. A more elaborate discussion of the spin arrange- 
ment will be given elsewhere. Here we only point out 
that (a) our calculations with the dimerized dimer lat- 
tice indicate that very large spin moments are possible 
when the intra-dimer hopping integrals are large com- 
pared to the inter-dimer hopping (see Fig. [To]), in quali- 
tative agreement with the observed very-large magnetic 
moment in ^BEDT-TTF) 2 Cu(CN) 2 CE3, and (b) each 
dimer of BEDT-TTF molecules has the cartoon occu- 
pancy of 10 or 01 and the ...1100... ordering along one di- 
rection and ...1010... ordering along another (see Fig. |^), 
thereby reducing the spin frustration among the dimer 
sites forming the triangular lattice. In the absence of 
this population difference within each dimer cell (and 
the population difference is a consequence only of dimer- 
ization of the dimer lattice) the frustration within the 
triangular lattice would have severely reduced magnetic 
moments. We further point out that a pseudo-gap in the 
spectrum of magnetic excitations has been observed in. 
the SDW phase of ^BEDT-TTJ) 2 Cu[N(CN) 2 ]BrE_T_ 
and k-(BEDT-TTF) 2 Cu(NCS) 2 _J; this is in agreement 
with the dimerization of the dimer lattice, since with- 
out the second dimerization there should be no spin gap 
within the 2D antiferromagnet. 

The material k-(BEDT-TTF) 2 Cu 2 (CN) 3 merits sepa- 
rate discussion. This material is not antiferromagnetic, 
and measurement of spin susceptibility due to the BEDT- 
TTF components exhibijtS| a steep drop below 10 K, sug- 
gesting SP-like behaviort_. This behavior is very similar 
to that in the BETS-based materials, which we discuss 
below, where we point out that for p = 1/2, this behavior 
is expected for the case of large t± (> tj_). 



E. A-(BETS) 2 GaBr z Cl 4 _„ (BETS = BEDT-TSF) 



These materials, discovered only recentlyEJ L_, are su- 
perconducting for < z < 0.8 and semiconducting for 
0.8 < z < 2.0. Thus the proximity between a semicon- 
ducting and a superconducting state that characterizes 
the TMTSF and the BEDT-TTF is also a characteris- 
tic feature of the A-BETS. In contrast to the TMTSF 
and the BEDT-TTF systems, however, the semiconduct- 
ing state-in the BETS is nonmagnetic and possesses a 
spin gap_J. Magnetic susceptibility studies indicate ab- 
sence of anisotropy in the susceptibility, and no spin-flop 
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transition (signature of antiferromagnetism) was found 
down to 10 K, which is close to the maximum supercon- 
ducting critical temperatuxe T c (onset 7.5 K, and even 
higher in certain samples)E3. The absence of the SDW 
is particularly perplexing here in view of the strong two- 
dimensionality predicted within extended Hiickel band 
calculations&l 

The lattice structures of the A-(BETS) 2 GaBr 2 Cl 4 _ z 
are knownEl The stacking of the organic donor molecules 
is very similar to the /3-BEDT-TTF systems, i.e., a 
nearly rectangular lattice with strong intrastack cou- 
pling, weaker transverse coupling, and very weak cou- 
pling along one diagonal. The nearly rectangular lattice 
permits comparison with our theory. One interesting fea- 
ture of the lattice structure is that the intrastack bonds 
have strengths that are W'SWS, exactly the structure 
expected for the r^ F ^ lattice in Fig. ||(b) as well 
as the dimerized dimer lattice. We believe that while 
the difference between the strong and weak bonds is a 
crystal structure effect, the further dimerization of the 
dimer lattice is a consequence of the BCDW instability 
discussed here. r— ■ 

Hartree-Fock calculations by Seo and FukuyamaEa 
within an anisotropic Hubbard Hamiltonian gave an an- 
tiferromagnetic ground state instead of the nonmag- 
netic state. Since Hartree-Fock calculations overesti- 
mate antiferromagnetism, these authors then chose the 
U — » oo limit of Hubbard model to arrive at a dimer- 
ized, anisotropic 2D Heisenberg spin Hamiltonian, each 
lattice site of which corresponds to one dimer of the 
original BETS lattice. The antiferromagnetic-SP bound- 
ary within the 2D dimerized Heisenberg spin Hamilto- 
nian has been investigated by Katoh and Imada using 
QMC simulations^. For the longitudinal and trans- 
verse exchange integrals derived by Seo and Fukuyama, 
the QMQjjalculations still predict the antiferromagnetic 
structural. Seo and Fukuyama explain the spin gap 
in A-BETS by claiming that the second dimerization of 
the dimer lattice (i.e., intermolecular distances W'SWS, 
instead of WSWS) takes these systems to the ID side 
of the 1D-2D antiferromagnetic-SP boundary, exactly as 
(TMTTF) 2 PF 6 , even though the actual transverse hop- 
ping integrals are large. 

We believe that the problem faced by these authors 
arises entirely from their effective 1/2-filled band approx- 
imation. As seen in Fig. |l^, the dimerization of the dimer 
lattice enhances the SDW in the region of small to in- 
termediate t± and therefore cannot be the origin of the 
spin gap or supposedly ID behavior. Recall also that 
(TMTTF) 2 PF 6 , which is certainly on the ID side of the 
1D-2D boundary, is nonsuperconducting. In contrast, A- 
BETS does become superconducting and that too at a T c 
that is considerably higher than that in the (TMTSF)^, 
indicating what we believe to be strongly 2D characters. 
We believe that the solution to this puzzle lies in recog- 
nizing the p =1/2 character of the (BETS) 2 X. An essen- 
tial difference between the effective 1/2-filled band model 
of Seo and Fukuyama and ours is that within the former, 



there are only two regions, nearly ID and 2D, with the 
spin states as singlet and antiferromagnetic, respectively. 
Our work indicates that there are three distinct regions, 
singlet, antiferromagnet, and singlet again, as a function 
of increasing t± , independent of whether one assumes a 
1/4- filled band or an effective 1/2-filled band. We there- 
fore believe that a more natural explanation of the spin 
gap phase is obtained within our theory, with the singlet 
ground state in semiconducting BETS not being due to 
t± that is too small, but due to a t± that is too large 
(> t c ) to give SDW. This would be in agreement, with 
the strong two-dimensionality of these systemsE-sE3. We 
believe that the same explanation also applies to the k- 
(BEDT-TTF) 2 Cu 2 (CN) 3 , discussed in the above. We 
predict that experiments that can probe charge ordering 
will find two kinds of BETS molecules with different elec- 
tronic populations, with greater charge densities on the 
two BETS molecules that are linked by the W bond. 



VII. POSSIBLE IMPLICATIONS FOR ORGANIC 
SUPERCONDUCTIVITY 

What might be the implications of our BCSDW and 
BCDW states to organic superconductivity, the mecha- 
nism for which remains unclear despite two decades of 
research? We present here several partial responses to 
this challenging question. 

First, given the the robustness of the BCDW/BCSDW 
in the exactly 1/4-filled band, we believe that the su- 
perconductivity must be the result of weak incommen- 
surability in the actual materials. Specifically, we sug- 
gest and discuss in more detail below, that supercon- 
ductivity arises from the pairing of commensurability 
defects in the background BCDW/BCSDW. That such 
weak incommensurability exists is strongly indicated by 

(i) the observation of a zero-energy mode in the optical 
conductivity^!^ of (TMTSF) 2 PF 6 and (TMTSF) 2 C10 4 ; 

(ii) the observatiorLxif Frohlich mode sliding transport in 
the same materialsEj; and (iii) the observation of a '^par- 
tially gapped Fermi surface" in the metallic regior£3 of 
a-(BEDT-TTF) 2 KHg(SCN) 4 . Extremely interesting ce- 
sults in this context were reported by Komatsu et. al.EJ, 
who showed that the superconductivity in k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 was due to a subtle change in the va- 
lence state of the Cu. The pure K-phase material is a 
semiconductor with the Cu valence of +1. According 
to the authors of Ref. ^3|, the superconducting phase 
corresponds to a different material («;' in the authors' 
notation) in which some of the Cu (several hundred 
ppm) have acquired valency 2+. This was confirmed 
from ESR studies. The increase in Cu valency decreases 
the overall negative charge on the anion, and there- 
fore the overall positive charge on the cation, provid- 
ing a weak incommensurability that appears to be essen- 
tial for superconductivity^. This result lends credence 
to our suggestion that organic superconductivity arises 
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from the pairing of commensurability defects within the 
BCDW/BCSDW background. 

Second, the similarities between the organic and high 
temperature oxide superconductors have been 
out in recent years by several research group; 
One obvious apparent similarity between these two 
classes of superconductors is the proximity of the SDW 
to superconductivity. Our studies suggest that supercon- 
ductivity in the organics is actually occurring at the in- 
terface of a Coulomb-induced BCDW that for a range of 
t± coexists with the SDW. It therefore seems more likely 
that the pair binding is actually driven by the BCDW, 
and not the SDW, although it is probable that the sym- 
metry of the pairing state may depend on the SDW (see 
below). As noted above, the experimental observation 
of superconductivity in the Ju(BETS)2GaCi4 (where no 
proximate SDW is observedEZI) supports this view. An 
important implication of this perspective is that it casts 
doubt on recent spin-fluctuation theories of organic su- 
perconductivitv-within the effective 1/2-filled correlated 
electron modeO~t2£|. The consequences of this conclu- 
sion from the organics for the high T c materials are un- 
clear, but it is perhaps not irrelevant in this context to 
point out that evidence for superconductivity within the 
2D nearly 1/2-filled Hubbard model, which for laxge-L/ 
has strongly AFM behavior, has remained elusi* 
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despite more than a decade of intense research!: 

Third, there are striking similarities between this 
"doped" BCSDW/BCDW scenario and several other the- 
oretical suggestions of superconductivity induced by dop- 
ing of exotic "paired" semiconductors. As we have noted 
previously, the BCSDW and the BCDW states are very 
similar to the "paired electron crystal" (as opposed to 
the monatomic Wigner crystal) found by Moulopoulos 
and Ashcroft for the intermediate density electron gasEZl. 
Superconductivity near the "melting" transition of the 
paired electron crystalJaaSpbejen conjectured by a number 
of authors in the pastJi23~li2a, even before the discovery 
of organic or high T c superconductivity. The commensu- 
rate BCDW is also qualitatively similar to a "negative U - 
positive V" effective 1/2-filled extended Hubbard model, 
with the effective lattice sites sites consisting of (a) the 
"occupied" pair ('1-1') of nearest neighbor sites, and (b) 
the "unoccupied" pair ('0-0') of nearest neighbor sites, 
in Fig. 0(c). Within this scenario, there is an effective 
attraction between the carriers on the "occupied" pair 
of dimer sites, but an effective repulsion between two 
pairs of occupied dimers. For models of this type, it 
is known that diagonal and off-diagonal long-range order 
can ii«wfip£iple coexist slightly away from commensurate 
fillingli3 _ IH3. Further, Imada has studiedti^l a 2D spin- 



Peierls state (not possible in the monatomic 1/2-filled 
band) in which each composite site is again a dimer, with 
the dimer sites now having occupancies '10' and '01' (see 
Fig. ||(b) and note that the bonds between a '10' and '01' 
and between a '01' and a '10' are different, giving rise to 
a spin-Pcicrls-likc behavior). His numerical simulations 
find evidence for superconductivity in the hypothetical 



doped 2D spin-Peierls statet!3. Finalls^-Emery, Fradkin, 
and Kivelson have recently suggestedEIJ that supercon- 
ductivity can exist for incommensurate fillings in models 
that support stripe phases and i n wh ich a spin gap is 
present. Since the analysis in Ref. 114 does not make di- 
rect contact with an initial microscopic Hamiltonian, but 
rather posits the form of the effective Hamiltonian in the 
vicinity of an unpinned stripe phase, it is not possible im- 
mediately to make detailed comparisons with our res ults. 
We can, however, make two comments. First, Rcf. 114 



reflects the widespread belief that models within which 
a spin gap persists in the doped state are strong candi- 
dates for a microscopic theory of correlated superconduc- 
tivity. Our preliminary numerical evidence suggests that 
both the BCDW and the BCSDW will continue to have 
a spin gap when doped; further work is in progress to 
confirm this. Second, regarding the attractive possibility 
that our BCSDW/BCDW state provides the background 
charge order within which commensurability defects may 
pair to form a superconducting state, we note that the 
occupancy schemes in Fig. 0(c) and Fig. ||(a) and (b) re- 
semble intersecting stripes, where each stripe is obtained 
by connecting the '1—1' bonds along the x- and x+y (- 
x+y) directions. 

The possible BCSDW/BCDW to superconductor tran- 
sition in the organic CTS clearly requires further study. 
We close our present discussion of this topic with com- 
ments on three important open issues: (i) the possible 
mechanism for superconducting pairing; (ii) the problem 
of phase separation; and (iii) the symmetry of the order 
parameter. 

First, the possible mechanism for pairing of commen- 
surability defects within the 2D BCDW can be visual- 
ized most simply in the rigid bond limit, where nearest- 
neighbor bonds retain their individual distortions inde- 
pendent of the occupancies of the sites linked by these 
bonds. The commensurate BCDW in this limit can 
be viewed as consisting of "quasimolecules" , where each 
quasimolecule is a "1-1" dimer. If two holes are now 
removed from the system, it is energetically preferable 
to destroy one "quasichemical bond," thereby creating 
an intersite (small) bipolaron, as opposed to destroying 
two bonds and creating two polarons. Thus, within the 
W'SWS structure (t s > tw > tw), each W bond acts 
as a "negative-U" center in the rigid bond limit. It is 
of course highly unlikely that superconductivity can be 
obtained, at least at the experimental T c , due to con- 



densation of small bipolaronstlj, so this might appear 



to present a serious problem for this proposed mecha- 
nism. In fact, when one goes beyond the oversimplified 
rigid bond limit to the full model that correctly reflects 
the cooperation between e-e and e-ph interactions in the 
1/4-filled band, one finds that the actual commensurabil- 
ity defects are more like the extended, "resonant" (and 
therefore mobile) bipolarons that are indeed candidates 
for explaining superconductivity in strongly correlated 
systemsE12lEi3. To understand this in detail, consider 
again the weakly incommensurate BCDW, starting from 
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the ID limit, but now with the e-ph interactions included. 
Below the 4kp transition temperature T^fi but above 
the 2k transition temperature T2k F , incommensurabil- 
ity leads to fractionally charged solitons wit h char ge e/2, 
and each vacancy creates two such defects! 1 15 H 116 I. Pre- 
vious work has assumed that the soliton charge remains 
e/2 even below the 2kp- transition, wiiich implies that 
two vacancies create four such defectstlEl. However, Ref. 
|116| assumes that the site charge density remains uniform 
even below the ID 2kp SP (dimerization of the dimer lat- 
tice) transition, which is precisely what we have shown 
here not to be the case. Indeed, as a consequence of this 
spatial charge inhomogeneity (charge ordering) , the "soli- 
tons" now acquire integer charge (i.e., two fractionally 
charged solitons bind to give a single solitop-yith charge 
+e), as we have shown explicitly elsewhereEITJ. A pair of 
added vacancies within the ID BCDW below T2fe F there- 
fore creates (only) two solitons. In the strictly ID limit, 
these do not bind, but with increasing t±, one expects 
binding to a large bipolaron. The source of this bind- 
ing is precisely the same as the source oLsoliton confine- 
ment in coupled chains of polyacetyleneEl: in the region 
between the two defect centers the phase relationships 
between the BCDW's on neighboring chains is different 
from the preferred one (viz., tt), and therefore a large 
separation between the defect centers would increase the 
energy (linearly with increasing separation) . There exists 
therefore a space- dependent interaction between the po- 
larons, which is repulsive at short range but attractive at 
some (i^-dependent) intermediate range. The bipolaron 
size, as well as its dimensionality, depends on t± (as well 
as on U and V). There is currently limited analysis of 2D 
large bipolarons in the strongly correlated limit, although 



some results suggest that these can indeed be mobildli^ 



Within this scenario, superconductivity occurs due to the 
condensation of these large bipolarons, w hich is not pre- 
cluded by the theoretical analysis of Ref. 111. Resolving 



the question of whether static distortion is sufficient, or 
whether dynamical phonons will have to be included, will 
require further work. 

Second, in many existing models of superconducting 
pairing involving correlated electrons, the interactions 
that bind two particles also lead to phase separation, 
since the attraction producing pairing does not satu- 
rate. Eerhaps the best known example of this is the t-J 
mode£2J : Ei3 away from 1/2-filling. In contrast, within 
any "negative U" model there does exist a saturation in 
this attraction (since a single site can at most have two 
electrons), and the analogy between our BCDW model 
and the effective 1/2-filled "negative U - positive V" case 
suggests that phase separation will also not occur here. 
Further, the immediately previous discussion of the pro- 
posed binding mechanism makes clear that with small 
but macroscopic (say, 1%) concentration of commensu- 
rability defects, there is no particular energetic advan- 
tage in creating additional polarons or bipolarons prox- 
imate to the original bipolaron (in contrast to, say, the 
t-J model, where there is such an energetic advantage). 



Third, what symmetry do we expect for the supercon- 
ducting order parameter in our model? This is clearly a 
challenging issue, particularly since even with the same 
BCDW background the pairing symmetry in the highly 
anisotropic TMTSF might be different from that in the 
more two-dimensional BEDT-TTF and BETS. Several 
recent experiments have presented evidence consistent 
with nodes p_tlie-.superconducting gap function in the 
BEDT-TTFE3TB. This is reminiscent of d-wave sym- 
metry of the superconducting order parameter in the 
high temperature copper oxide based superconductors. 



On the oth 



t-hand, Lee at al. have recently presented 



pviHpnrpl 123 H 124 l suggesting that a spin triplet p-wave pair- 
ing is necessary to explain data in (TMTSF^PFg, where 
the upper critical field H C 2 shows no saturation with 
the field in the plane of the organic molecules and ex- 
ceeds the Pauli paramagnetic (Clogstpn) limit expected 
to hold for singlet superconductoralH3 and the tempera- 
ture dependent Knight shift measurements of 77 Se show 
that the spin susceptibility remains unaltered through 
the sup erco nducting T c li£fl. Within the continuum RG 
theorieal^O triplet superconductivity does indeed occur 
proximate to the SDW. However, within the discrete ex- 
tended Hubbard model, triplet superconductivity occurs 
within a very narrow region of the the "positive U - neg- 
ative V" sector of the U — V phase diagram, bounded by 
the SDW phase and a phase segregated phased. Triplet 
pairing thus will not only require a change in sign of the 
nearest neighbor Coulomb interaction within our origi- 
nal Hamiltonian of Eq. (1), but will also occur for a very 
narrow critical range of this parameter. But to resolve 
definitively the issue of the symmetry of the order pa- 
rameter within our model will be a non-trivial task, as 
the consequences of the interplay between e-e and e-ph 
interactions, as well as the effects of anisotropy, must be 
properly understood. 
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APPENDIX 1: THE AFM-SINGLET 
TRANSITION FOR WEAK ANISOTROPY 

Our goal here is to understand the second AFM-to- 
singlet transition that should occur in the quarter-filled 
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band for large t± from a perspective that is different from 
the one presented in section III. Specifically, we refer to 
the antiferromagnetic dimer lattice of Fig. |l](b) with weak 
intrachain interdimcr links, and the frozen valence bond 
state of Fig. 0(c), in which one of the interdimer links 
(W in the notation of section III) is now stronger than 
the other (W in the notation of section III), and is a 
singlet bond. We aim to give variational arguments at 
the simplest level that (a) point out the difference be- 
tween p — 1 and p — 1/2, and (b) indicate that the 
frozen valence bond state of Fig. pYc) dominates over 
the antiferromagnetic dimer lattice of Fig. |l|(b) for large 
ij_ and therefore the dimerization of the dimer lattice is 
unconditional. The argument given below is not to be 
considered as a proof, but rather, it provides convincing 
physical motivation for the numerical work discussed in 
section IV. 

Note that our discussion here is limited to the relative 
stabilities of two insulating states, and not the competi- 
tion with any metallic state. We consider only the simple 
Hubbard Hamiltonian with V = (since for p — 1 /2 the 
periodicity of the CDW is the same for all V < V c and 
while for p — 1 the V merely reduces the effective on-site 
correlation) for t± = 1. For completeness we begin by 
repeating the variational argument for the dominance of 
the SDW over the BOW in p = 1. Consider the Heisen- 
berg antiferromagnetic spin Hamiltonian 

H=jJ2Si-S j (11) 

m 

Consider also the singlet variational state (1,2)(3,4)....(N- 
1,N), with singlet bonds between nearest neighbors in ID 
and the Neel state ... Till The energy of an isolated 
singlet bond is -(3/4) J while that of a 2-site Neel state is 
-(1/4) J. The overall variational energy of the singlet state 
in ID is -(3/8)NJ and that of the Neel state -(1/4)NJ, 
so that the singlet dominates over the Neel state in ID. 
In the 2D isotropic N x N lattice, we compare (1) the 
frozen valence bond state in which each chain still has 
the same spin couplings as in ID (note that at the level 
of our approximation the relative phases between con- 
secutive chains make no difference), and (2) the 2D Neel 
state. The variational energy of the frozen valence bond 
state is-(3/8)N 2 J, but now because of the larger number 
of nearest neighbors the energy of the Neel state has a 
lower value -(1/2)N 2 J, which therefore dominates over 
the frozen valence bond state. Thus for p = 1 in ID the 
SP state dominates, while in 2D the SDW wins over the 
SP state. While this argument may appear simplistic, it 
nevertheless predicts the dominance of the antiferromag- 
net over the singlet in p = 1. 

Consider now the isotropic 2D dimerized p = 1/2 lat- 
tice with moderately strong dimerization (Fig. 1(b)). 
The effective 1/2-filled band is clearly a SDW, with 
the dimerization pattern being necessarily "in-phase" be- 
tween consecutive chains, as shown in Fig. [j](b) (to pre- 
vent confusion in what follows we have not shown the 



bonds in Fig. |l|(b), but a strong bond between the two 
sites within the parentheses and weaker interdimer bonds 
have been assumed) The individual site populations are 
equal in this state and each is exactly 1/2. Our con- 
tention is that this state has a higher variational energy 
than that reached by further dimerization of the dimer 
lattice, which gives the p =1/2 frozen valence bond state 
shown in Fig. |l|(c), where there occur interdimer singlet 
bonds and site occupancies ...1100... (the singlet bonds 
in Fig. 0(c) are between the "occupied" sites) . The rea- 
son for this is that unlike in p — 1, the exchange inte- 
grals that describe the effective Heisenberg models in the 
SDW and the singlet are now different, in spite of the fact 
that both Heisenberg systems are derived from the same 
Hubbard Hamiltonian. In Fig. 0(c), we are considering 
isolated singlet bonds, with site occupancies of 1. and J 
is clearly 2t 2 /U, exactly as for p = 1. In Fig. 0(b), on 
the other hand, the exchange integral has to correspond 
to a true p = 1/2 system, since each site occupancy is 
now 1/2. The exchange integraLJ' for arbitrary p in 
ID is 2{t 2 /U)p[l - sin(2Trp/)2npf3, so that for p= 1/2 
we have J' = (1/2) J along each chain (the x-direction). 
This expression is strictly true only in the ID undistorted 
chain, and for the distorted ID chain or in 2D one needs 
to calculate J' from comparing singlet-triplet gaps within 
the structure corresponding to Fig. 0(b) and within the 
1/2-fillcd band. We have calculated these gaps for fi- 
nite lattices separately for the longitudinal and trans- 
verse directions and have found that while J' = (1/2)J 
is quite accurate for the longitudinal direction, the J' in 
the transverse direction is even smaller (the difference 
between the longitudinal and transverse directions origi- 
nates from dimerization along the longitudinal direction 
only), with the restriction that only interdimer hops lead 
to spin exchange. Even if we consider the largest possible 
value for J' = J/2, the variational energy of the Neel state 
in Fig. 0(b) is then -(1/2)(N 2 /2)J' = -(1/8)N 2 J, while 
that of the frozen valence bond state in Fig. 0(c) (with 
N 2 /4 singlet bonds) is -(3/4)N 2 /4J = -(3/16)N 2 J. Thus 
the frozen valence bond state dominates over the dimer 
SDW, implying that the dimerization of the dimer lat- 
tice is unconditional, and the difference from the simpler 
p = 1 case arises from the smaller (by factor of 2) ex- 
change integral in the uniform dimer lattice of Fig. 0(b). 
The above approach is obviously simplistic, but no more 
so than the physical argument for the dominance of the 
SDW in p=l. 

APPENDIX 2: NUMERICAL METHODOLOGY 

A. Contstrained Path Monte Carlo (CPMC) 

The j-CPMC ground-state quantum Monte Carlo 
methocS uses a constraining trial wavefunction to elim- 
inate exponential loss of signal due to the Fermion sign 
problem. Although the method has been thoroughly 
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bench-marked against known results for the Hubbard 
model, the method is non-variational, and it is impor- 
tant to check its accuracy in every new system against 
exact results and to use a variety of different trial wave- 
functions. The current application is different from pre- 
vious ones in including the V interaction, as well as 
in the choice of the bandfilling (previously tested cases 
were for bandfillings close to 1/2). Furthermore, previ- 
ous work has shown that most accurate results are ob- 
tained when the trial noninteracting wavefunctions have 
"closed-shell" nondegenerate configurations. In the next 
subsection it is shown that the proper boundary condi- 
tions for simulating coupled 1 / 4-fillcd band chains in- 
volves having 4n electrons (where n is an integer) per 
chain. This implies degeneracy of the trial wavefunctions 
at t± = and again near t± = 1. It is thus necessary to 
check the accuracy of the method for our purpose, and 
this was done by comparing CPMC results with the exact 
results for the 8x2 lattice. 

Fig. [l3] summarizes the results of the bench-mark en- 
ergy calculations for an 8x2 lattice, periodic in the x- 
direction, with U — 6 and V = 1. Both undistorted 
and and the 2kp (r 2 k F ^ 0,r4fc F = in Eq. (7)) bond 
distorted systems were compared, where for the the uni- 
form lattice all hopping integrals were taken to be 1.0, 
while for the distorted system they were 1.14, 1.0, 0.86 
and 1.0 (as in Fig. 2(a)). For this amplitude of the 2ki? 
distortion, 
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FIG. 13. Percent errors in the CPMC energies for (a) 
undistorted and (b) 2k^ bond-distorted (the hopping inte- 
grals in the distorted lattice correspond to those in Fig. ^|(a)) 
8x2 lattices with (7 = 6 and V = 1. Triangles are for the 
free-electron trial function; circles for the UHF trial function. 



the absolute value of AE is only 0.3% of the total 
energy (at t± = 0.4). Such a small energy difference 
is not easy to measure within quantum Monte Carlo. 
We note that energy differences of this order of mag- 
nitude have also been calculated using CPMC to. study 
hole binding in the the 3-band Hubbard modelEa. The 
CPMC values are scaled for At -> from At = 0.05 
and At = 0.1 to remove the Trotter discretization er- 
ror. The trial wavefunctions used were either the free- 
electron wavefunction, or an Unrestricted Hartree-Fock 
(UHF) wavefunction with U = 2 and V = 0.5. Hartree- 
Fock wavefunctions with larger U and V gave less ac- 
curate results, probably due to the tendency of UHF to 
exaggerate AFM correlations. In Fig. [l3| the UHF trial 
functions produced larger errors than the free-electron 
trial functions for the distorted system at small t± be- 
cause the SDW correlations there are exaggerated by the 
UHF approximation. The CPMC systematic errors are 
largest at small t±_ (< 0.2) and large t± (> 0.8) possibly 
due to the degeneracies in the one-electron occupancies 
at t± = and t± = t. However, at large t±, the UHF 
trial wavefunction produced slightly more accurate re- 
sults for the 8x2 distorted lattice possibly because the 
numerically-derived UHF wavefunction breaks some of 
the symmetry of the non-interacting wavefunction. In 
the intermediate t± ~ 0.4 regime, the CPMC energies 
are indistinguishable from the exact energies within the 
statistical error. The accuracy of the CPMC method in 
this region is very reassuring, since for the noninteracting 
case, at t± = 0.4 the distortion has already vanished. 

In addition to comparing energies, we have also com- 
pared charge densities and spin-spin correlation func- 
tions. Table I compares the charge densities and spin- 
spin correlations computed by CPMC for the 8x2 dis- 
torted lattice at t± = 0.4. The agreement with the exact 
result is not as good as for the energy (typically 1-5% for 
the charges and 5-10% for the spin-spin correlations), but 
is more than adequate to identify the presence and peri- 
odicity of the broken symmetry states. Thus in general, 
we find the CPMC results are close to the exact results 



(p 3 ) 


(44) 


j 


exact 


CPMC 


U 


exact 


CPMC 


1 


0.4799 


0.4756(6) 


1,9 


-0.06095 


-0.0585(7) 


2 


0.5201 


0.5250(6) 


1,10 


-0.03215 


-0.0312(7) 


3 


0.5201 


0.5240(6) 


1,11 


0.01408 


0.0161(7) 


4 


0.4799 


0.4772(6) 


1,12 


-0.02698 


-0.0231(6) 








1,13 


-0.07299 


-0.0687(6) 








1,14 


-0.03085 


-0.0268(5) 








1,15 


0.01408 


0.0158(6) 








1,16 


-0.02552 


-0.0239(7) 



Table I. Comparison of CPMC and exact charge density and 
spin-spin correlations for an 8x2 system with U = 6, V = 1, 
tj_=0.4, with the same distortion of hopping integrals as in 
Figure (|l3|). Sites on the first chain are numbered 1-8, those 
on the second chain 9 - 16. 
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for both energies and correlation functions, except for 
very small or large t±. 



B. Boundary Conditions 

As noted above, we determine the proper combina- 
tions of lattices and boundary conditions for the nu- 
merical simulations by the requirement that nonzero t± 
destabilizes the BCDW for noninteracting electrons with 
those boundary conditions on that particular finite lat- 
tice: i.e., we require the finite lattices to reflect correctly 
the known behavior of the noninteracting case in the ther- 
modynamic limit. 

Consider an N x M lattice, with N sites along the 
chain and M chains. To avoid odd/even effects, consider 
an even number of electrons per chain. This number can 
then be either 4n or 4n + 2, where n is an integer. To 
obtain a f/4-filled band, one can then have N = 4n x 
2 = 8n or N = (4n + 2) x 2= 8n + 4. The proper N 
for our purpose is N = 8n (i.e., 4n electrons per chain). 
This follows from the one-electron energy levels of cou- 
pled chains with 4n and 4n + 2 electrons per chain. In 
Fig. |l4] below we have shown the one-electron energy lev- 
els for the undistorted 8x2 (top panel, labeled a) and 
12 x 2 (bottom panel, labeled b) lattices (both periodic 
in the x-direction), corresponding to t± = on the left 
and 0.1 on the right in both cases. In the 8x2 lat- 
tice, the degeneracy at t± = will lead to spontaneous 
distortion. For nonzero t± and a 7r-phase shift between 
chains (which gives lower energy than phase shifts of 
or S), the pairs of one- electron levels that are coupled by 
phonons with wave- vector (2kf, tt) are (- 0) and (+ 
2£, tt); and (+ ?f, 0) and (- ?f, tt). The finite gap that 
occurs for t± ^ between each pair of one-electron 
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FIG. 14. Occupancies of the one-electron levels for the 
undistorted (a) 8 x 2 lattice, with t±_ — (left) and t± = 
0.1 (right) and (b) 12 x 2 lattice, also with t± = (left) and 
tj_ = 0.1 (right). 
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FIG. 15. One-electron AE/AE versus t± for (a) N = 8n 
and (b) N = 8n+4. In each case results for several N x M 
lattices are shown. 

levels coupled by the (2kf, tt) phonon indicates ab- 
sence of nesting and the destabilization of the distortion. 
This energy gap increases with t±, leading to a decrease 
in AE with t± for N = 8n (see Fig. 15 for details), as 
occurs in the thermodynamic limit. In contrast, consider 
the 12 x 2 lattice, in which the one-electron ground state 
is non-degenerate. There is now a nonzero energy gap 
between the levels coupled by the 2kp electron-phonon 
interaction already at t± = (k x — - ?5 and k x = + 
k x = + and k x — - te). With nonzero tj_, and 
once again a 7r-phase shift between the chains, the en- 
ergy gap between the levels (- y|, tt) and (+ y|, 0), and 
similarly that between the levels (+ ff, tt) and (- 0), 
decreases, indicating that the tendency to distort here 
increases with inter-chain coupling, at least for small to 
moderate t±. 

For large N, the difference between N = 8n and N = 
8n+4 vanishes, as is shown in Fig. where Figs, [l^ (a) 
and (b) show the behavior of AE(t±) for N = 8n and 
8n+4, respectively. The qualitative behavior (destabi- 
lization of the distortion) is the same for all N = 8n, and 
monotonically decreasing AE is also seen for N = 8n+4 
for large N, but finite size effects (increasing AE at small 




FIG. 16. One-electron AE/AE for the 8 x 4, 8 x 8, 16 x 
4, 16 x 8 and 16 x 16 lattices for the 2kF bond distortions 
as in Fig. 0(a). 
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to intermediate t±) are strong even for N = 28, a chain 
length already too large for accurate 2D many-body cal- 
culations. The correct qualitative behavior of all N = 8n 
is the basis of our choice of these N. 

In contrast to the choice of N, there is no immediate 
restriction on the choice of M, the number of chains, ex- 
cept that M should be even, to avoid even/odd effects. 
M = 4n and 4n + 2 both show the same qualitative be- 
havior, as seen from the plots of AE versus t± in Fig. [If], 
for several M = 4n lattices (M = 4n + 2 are included in 
Fig. |l5|) . Thus both M = 4n and 4n + 2 are appropriate. 
Our choice of M = 4n + 2 is based on two reasons. First, 
exact diagonalization calculations on the 8x2 lattice al- 
lows comparisons to results obtained within CPMC, and 
the exact diagonalizations cannot be done for the next 
larger appropriate lattice, viz. 8x4. Second, the M = 
4n lattices are characterized by one-electron Fermi level 
degeneracies for t± ^ (even though the degenerate lev- 
els are not coupled by (2kp, ir) phonons), and the absence 
of a single well-defined one-electron wave-function would 
make the CPMC calculations considerably more difficult 
than for M = 4n + 2 lattices, which have non-degenerate 
one-electron levels for nonzero t±. 



C. UHF calculations of bond distortion 



values of U and V gave better results when used as 
a CPMC trial function (compared to a numerically ex- 
actly solved 8x2 system). Fig. [l?] shows the normalized 
energy gain from a 2k^ distortion for two different lat- 
tices, within the UHF approximation. The UHF results 
show that AE/AEq remains close to 1 for at least up to 
tj_ ~ 0.4, indicating a tendency to persistent distortion 
up to this t±. Although AE/AEq begins to decrease at 
still larger t ± , these calculations are for a relatively small 
value of U, and as discussed in section III, the range of t± 
over which the distortion should persist increases with U. 
Thus the qualitative effects of the e-e interaction are al- 
ready visible within the UHF approach at small U, while 
a fully persistent broken symmetry state will occur only 
for larger values of the e-e interaction that are beyond the 
scope of the UHF. Given that the UHF approximation 
predicts a vanishing of the bond dimerization in the 1/2- 
filled band for a fairly small U c (the actual magnitude 
of U c depends on a), in contrast to the correct result 
that there is an enhancement of the dimerizationlil for 
< U < 4, the present results, showing a persistence 
of the distortion for moderate t±, is initially perplexing. 
The reason for the correct prediction in this case is that 
the UHF exaggerates the SDW, which destroys the BOW 
in the 1/2-filled band, but has a co-operative interaction 
with the 1 / 4-filled band BOW for small to moderate t± . 



As discussed in the subsection on methods in this Ap- 
pendix, UHF trial wavefunctions for the CPMC calcu- 
lations were constructed for regions where one-electron 
wavefunctions were degenerate. Since UHF calculations 
give reasonably correct results in the small U, V range it 
is also of interest to determine the tendency of the 2D 
lattice to distort within the UHF approximation. One 
advantage of this procedure is that much larger lattices 
than those discussed in section IV can be tested. We re- 
port these results here. We have chosen relatively small 
U and V for two reasons: the UHF procedure does not 
converge well for larger interactions, and the smaller 
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FIG. 17. AE/AEo versus t± for a 2kf bond distortion 
( r 4k F = 0) for noninteracting and the interacting lattices 
within the UHF approximation. Intra-chain hopping integrals 
for the distorted lattices are as indicated in Fig. 0(a). 



APPENDIX 3: SPIN CHARACTER OF THE 
GROUND STATE 

As discussed in Appendix 2, the proper boundary con- 
dition for the numerical evaluation of the electronic en- 
ergy gained upon bond or site distortion in p = 1/2 in- 
volves finite N x M lattices with N = 8n. This requires 
the number of electrons per chain to be 4n, and it is 
known that in ID periodic undistorted rings with p ^ 1, 
the ground state has overall spin S = 1 instead of for 
any nonzero U. 



t± 


undistorted 


2kp distortion 




S=0 


S=l 


S=0 


S=l 


0.01 

0.025 

0.05 


-9.335651 
-9.337570 
-9.344122 


-9.335637 
-9.336944 
-9.341546 


-9.352522 
-9.354380 
-9.361083 


-9.352228 
-9.353739 
-9.358425 



Table II: The S=0 and S=l energies of the 8x2 undistorted 
and 2kF bond-distorted lattice for U — 6 and V = 1. The 
lowest energy is S=0 for both undistorted and distorted cases. 

The spin of the ground state of the distorted periodic 
ring depends on its size and the magnitude of the Hub- 
bard U. For the values of the correlation parameters and 
bond distortion parameter in Fig. ^, the ground state in 
the N = 8 distorted periodic ring has S = 1, while the N 
= 16 ground state has S = 0. Thus the AE in Fig. | 
for nonzero e-e interaction corresponds to AEtt (i.e., 
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the energy gained by the triplet state upon bond distor- 
tion) for N = 8, and to AEts (undistorted state in S 
= 1, distorted state in S = 0) for N = 16. Whether or 
not the comparisons of the zero and nonzero ij_ are then 
meaningful is an important question. We present here 
the detailed results of three different sets of calculations, 
each of which indicates that our interpretation of the re- 
sults of Fig. |] (viz., strong tendency of the interacting 
1/4-filled lattice to distort at arbitrary t±) is correct. 

First, we have calculated the exact ground states of the 
8x2 lattice for t± as small as 0.01. In Table II we have 
given the S = and S = 1 energies of the 8x2 lattice for 
[7 = 6 and V = 1, for three small values of t± . The cou- 
pled chain system is in the S = state for both zero and 
nonzero bond distortion for the smallest nonzero t±_ . The 
important point now is that instead of choosing the single 
isolated chain as the standard in Fig. [|, we could have 
also chosen the coupled chain system with tj_ = 0.01 as 
the standard, provided the distortion of the t± = 0.01 lat- 
tice is also unconditional. Even if the nesting ideas were 
valid, we believe that the coupled chain system with t± 
= 0.01 is unconditionally distorted and then the results 
in Table II clearly show that AE increases with further 
increase in t±, indicating enhanced distortion relative to 
tj_ = 0.01. The error bars in the CPMC calculations pre- 
vent us from performing similar calculations for the 8 x 
6 or the 16 x 6 lattices, but the overall similarities in the 
(a) occupancies of the one-electron levels for nonzero t± 
and (b) AE behavior, especially in the region t± < 0.4, 
preclude different behavior at small nonzero t±. 

We performed a second set of calculations for the 8 x 
2 lattice for very small values of U (with V = 0). Note 
that if the persistent distortion implied in Fig. ^ were 
merely due to our choosing the wrong reference point t± 
= 




FIG. 18. (a) AE vs. t± for the 8x2 lattice at small U 
(normalized to the value at t± — 0.05). Note the decrease in 
the AE. (b)AE/AE vs. tj_ for the 8x2 lattice at U = 100. 



(since exactly at this point AEq — AEtt), an appar- 
ently enhanced distortion for nonzero t± should occur 
for all nonzero U (since the single chain is S = 1 for all 
nonzero [/, while the coupled chain system has S = for 
all nonzero t± and U). On the other hand, if the results 
in Fig. |J are due to the confinement effect discussed in 
section III.D, then enhanced/persistent distortion should 
occur only above a threshold e-e interaction: for weak e-e 
interaction the behavior should resemble that of the non- 
interacting lattice (with enhanced or persistent distortion 
occurring for a small range of t± near t± =0). We show 
here the results of calculations at small U for site distor- 
tion (as opposed to bond distortion), since we also report 
calculations for very large U below, and the bond distor- 
tion pattern (the magnitude of r^p) is [/-dependent, but 
the site distortion pattern is not. The distorted lattice 
here has site energies + e, + e, - e, - e (with e = 0.1) 
over four consecutive sites, and a 7r-phase shift between 
the two periodic rings. Since the 2k^ CDW has a syner- 
getic coexistence with both the r4fc F ps BOW (Fig. |2|(a)) 
and the r 4fcF ^ BOW (Fig. |(b))Ell a persistent CDW 
also implies persistent BOW; we have confirmed this by 
calculating the expectation values of the bond orders. In 
Fig. 18(a) we show the AE behavior as a function of t± 
for both U = 0.5 and U = 1. Decreasing AE with t± 
is a clear signature that the tendency to distortion here 
decreases with increasing two dimensionality, since con- 
finement at these small U is not sufficient to give persis- 
tent distortion. Even though these calculations are with 
fixed site energies, the expectation values of the charge 
densities depend on t±, and our calculated CDW am- 
plitudes decrease with t±, as expected from Fig. |18|(a). 
This behavior is exactly opposite to that in Fig.lfWb), 
indicating again a decrease in distortion with t± at small 
U . Finally, we emphasize that similar calculations have 
also been done with fixed 2k^ bond distortion, and once 
again we observe decreasing AE and CDW amplitude 
with increasing t±. 

We performed a third set of calculations with very large 
U = 100, again with the same site distorted lattice but 
now with e — 0.2, since at this very large U, the energy 
gained upon distortion for e = 0.1 is very small. The re- 
sultant BOW here has strong 4kp component {r^p 0), 
and this is why the distorted lattice was chosen to be the 
2kp CDW in this and the above calculations, such that 
meaningful comparisons between these extreme cases can 
be made. At this large U, the energy difference between 
S = and S = 1 states is negligible. For example, for the 
ID 8-site periodic ring AE$s (electronic energy gained 
in the S = subspace, with both undistorted and dis- 
torted states in S = 0) = 0.06222, while AE TT (elec- 
tronic energy gained in the S — 1 subspace, with both 
undistorted and distorted states in S = 1) = 0.06224. 
Fig. [lj(b) shows the AE behavior as a function of t± 
(with e = 0.2 now). An enhanced CDW (and therefore 
BOW) is seen from as a function of t±, where the sin- 
glet and triplet data points at t± = are the same. As 
seen in Fig. ft8](b), the AE for nonzero t± is weakly en- 
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hanced now even when compared to AEss at t± =0. 
Once again, the behavior of the CDW amplitude is in 
complete agreement with the prediction from Fig. |l^(b) , 
viz., a weak enhancement of the CDW amplitude with 
t±. 

Considering the above three different sets of results, 
we therefore conclude that the results in Fig. ||, Fig. ||(b) 
and Fig. ^ are not artifacts, and the persistent distortion 
is real and a true confinement effect, as would also be 
expected from the "variational" arguments in Appendix 
1. 



D. Baeriswyl, D.K. Campbell and S. Mazumdar in Conju- 
gated Conducting Polymers, edited by H. Kiess (Springer 
Verlag, Berlin, 1992) pp. 7 - 133. 

2 A.J. Heeger, S. Kivelson, R.J. Schrieffer and W.-P. Su, 
Rev. Mod. Phys., 60, 781 (1988). 

3 S.N. Dixit and S. Mazumdar, Phys. Rev. B29, 1824 
(1984). S. Mazumdar and D.K. Campbell, Phys. Rev. 
Lett. 55, 2067 (1985). 

4 S. Kivelson, Phys. Rev. B 28, 2653 (1983). 

5 W.-P. Su, Solid St. Commun. 48, 479 (1983). 

6 J.E. Hirsch, Phys. Rev. Lett. 53, 2327 (1984). 

7 E. Manousakis, Rev. Mod. Phys. 63, 1 (1991) and refer- 
ences therein. 

8 S. Mazumdar Phys. Rev. B 36, 7190 (1987) and ibid, 39, 
12324 (1989). 

9 F.C. Zhang and P. Prelovsek, Phys. Rev. B 37, 1569 
(1988). 

10 S. Tang and J.E. Hirsch, Phys. Rev. B 37, 9546 (1988) 
and ibid, 12327 (1989). 

11 A.W. Sandvik, Phys. Rev. Lett. 83, 3069 (1999). 

12 I. Affleck and B. Halperin, J. Phys. A 29, 2627 (1996). 

13 Z. Wang, Phys. Rev. Lett. 78, 126 (1997). 

14 E. Dagotto and T.M. Rice, Science, 271, 618 (1996) and 
references therein. 

15 F.D.M. Haldane, Phys. Lett. 93A, 464 (1983). 

16 J. Solyom, Adv. Phys. 28, 201 (1979). 

17 V.J. Emery in Highly Conducting One- Dimensional 
Solids, edited by J.T. Devreese et al., (Plenum, N.Y. 
1979). 

18 J. Voit, Rep. Progr. Phys. 58, 977 (1995). 

19 J.E. Hirsch and D.J. Scalapino, Phys. Rev. B 27, 7169 
(1983) and ibid 29, 5554 (1984). 

20 S. Mazumdar, S.N. Dixit and A.N. Bloch, Phys. Rev. B 
30, 4842 (1984). 

21 We exclude from our discussions the extensive work on 
doped 2- and 4-leg ladder systems. It is not clear currently 
whether or not these results are relevant for the true 2D 
lattice. 

22 A true metallic phase may not be obtained because of dis- 
order effects, but this is a different issue that is unrelated 
to the question of the existence of density waves, and is 
hence ignored here. 

23 H.H. Lin, L. Balents and M.P.A. Fisher, Phys. Rev. B 56, 



6569 (1997). 

S. Mazumdar, S. Ramasesha, R. Torsten Clay and D.K. 
Campbell, Phys. Rev. Lett. 82, 1522 (1999). 
An abbreviated discussion of the evolution of the BCSDW 
into the BCDW for large t± can b e found in S. Mazum dar, 
R. T. Clay, and D. K. Campbell, |cond-mat/9910164 . 
K.C. Ung, S. Mazumdar and D.K. Campbell, Solid St. 
Commun. 85, 917 (1993). 

K.C. Ung, S. Mazumdar and D. Toussaint, Phys. Rev. 
Lett. 73, 2603 (1994). 

J. Riera and D. Poilblanc, Phys. Rev. B 59, 2667 (1999). 
T. Holstein, Ann. Phys. 8, 325 (1959). 
D. Jerome and H. Schulz, Adv. in Phys. 31, 299 (1982). 
T. Ishiguro, K. Yamaji and G. Saito, Organic Supercon- 
ductors, 2nd Ed., Springer Series in Solid State Sciences, 
Vol. 88 (Springer Verlag, Berlin, 1998). 
D. Jerome, Science 252, 1509 (1991). 
P.W. Anderson, Phys. Rev. Lett. 67, 3844 (1991). 
D.G. Clarke, S.P. Strong and P.W. Anderson, Phys. Rev. 
Lett. 73, 1007 (1994). 

Given that the issue of confinement in interacting electron 
systems has been hotly debated, it is pertinent to empha- 
size that we use this expression in a considerably narrower 
sense than has sometimes been used in the literature. Our 
interest lies strictly in the interacting 1/4- filled band case, 
and it is not being implied that the effect discussed here 
applies to other bandfillings. 

J.P. Pouget and S. Ravy, J. Physique I, 6, 1501 (1996). 
C. Bourbonnais, Synth. Metals 84, 19 (1997). 
K. Kanoda, Hyperfine Interactions, 104, 235 (1997). 
H. Kino and H. Fukuyama, J. Phys. Soc. Jpn. 65, 
(1996). 

R.H. McKenzie, 
(1998). 

P.M. Chaikin, J. 
L.P. Gor'kov, J. 



2158 



Comments Cond. Mat. Phys. 18, 309 



Phys. I. (France) 6, 1875 (1996). 
Phys. I. (France) 6, 1697 (1996). 

43 In all cases X are inorganic anions that play only indirect 
roles in accepting electrons from the organic cations, and 
do not have any active role in charge transport, although 
they do influence the magnitudes of the interchain hopping 
integrals. 

44 J.P. Pouget and S. Ravy, Synth. Metals, 85, 1523 (1997). 

45 S. Kagoshima, Y. Saso, M. Maesato, R. Kondo and T. 
Hasegawa, Solid St. Commun. 110, 479 (1999). 

46 The BCSDW has also been found in ID Hartree-Fock 
calculations within the dimerized dimer lattice, see N. 
Kobayashi and M. Ogata, J. Phys. Soc. Jpn. 66 3356 
(1997). The SDW component here, however, is an arti- 
fact of the Hartree-Fock approximation as in ID the spin 
couplings are necessarily singlet. 

47 The underlying x <-> y symmetry in the isotropic 2D 
limit implies that there are two degenerate orthogonal 2D 
BCDW states for t ± = t. 

48 R. T. Clay, A. W. Sandvik, and D. K. Campbell, Phys. 
Rev. B 59, 46 65 (1999). 
M. Nakamura, 



VI 



cond-mat/9909277 



50 J. Hubbard, Phys. Rev. B 17, 494 (1978). 

51 Note that Hartree-Fock calculations will predict this 
...1010... ordering for any non-zero V, and indeed will pre- 
dict the n-dimensional generalization of this order for an n- 



27 



dimensional hypercubic lattice. This artifact of the mean- 
field approach is known rigorously to be incorrect in ID 
and is shown by our later results to be incorrect in 2D. 
Only in the limit U » V » t does this "checkerboard" 
order actually become the dominant broken symmetry. 

52 D.J. Klein and W.A. Seitz, Phys. Rev. B 10, 3217 (1974). 

53 P. Jordan and E. Wigner, Z. Phys. 47, 631 (1928). 

54 J.D. Johnson and B. McCoy, Phys. Rev. A 6, 1613 (1972). 

55 H. Q. Lin, E. R. Gagliano, D. K. Campbell, E. H. Fradkin, 
and J. E. Gubernatis, in The Hubbard Model, edited by D. 
Baeriswyl et al. (Plenum Press, New York, 1995), pp. 316 
- 336. 

56 F.D.M. Haldane, Phys. Rev. Lett. 45, 1358 (1980); J. 
Phys. C 14, 2585 (1981). 

57 K. Moulopoulos and N. Ashcroft, Phys. Rev. B 48 11646 
(1993). 

58 S. Zhang, J. Carlson and J.E. Gubernatis, Phys. Rev. 
Lett. 74, 3652 (1995); Phys. Rev. B 55, 7464 (1997). 

59 M. Guerrero, J.E. Gubernatis and S. Zhang, Phys. Rev. 
B 57, 11980 (1998). 

60 S. Mazumdar, D. Campbell, R.T. Clay and S. Ramasesha, 
Phys. Rev. Lett. 82, 2411 (1999). 

61 K. Yamaji, J. Phys. Soc. Jpn. 55, 860 (1986). 

62 B. Domoulin, C. Bourbonnais, S. Ravy, J. P. Pouget and 
C. Coulon, Phys. Rev. Lett. 76, 1360 (1996). 

63 D.S. Chow, F. Zamborzsky, B. Alavi, D.J. Tantillo, A. 



Baur, CA. Merlic and S.E. Brown, |cond-mat/0004106 . 

64 H.J. Schulz in Low Dimensional Conductors and Super- 
conductors, edited by D. Jerome and L.G. Caron (Plenum, 
N.Y. 1987), pp. 95 - 112. 

65 C. Coulon, P. Delhaes, S. Flandrois, R. Lagnier, E. Bon- 
jour, J.M. Fabre, J. Physique 43, 1059 (1982). 

66 S.S.P. Parkin, J.C. Scott, J.B. Torrance and E.M. Engler, 
J. Physique Colloque, 44, C-3, 1111 (1983). 

67 In order to explain the "purely electronic CDW" , 
Kobayashi et al., in a ID calculation have suggested 
that the transition to the ...1100... CDW is driven by a 
second neighbor intra-chain Coulomb interaction V2 [N. 
Kobayashi et al., J. Phys. Soc. Jpn. 67, 1098 (1998)]. The 
synergic coexistence of the BOW and the CDW found by 
us indicates that a BOW must exist even in this case. 

68 G. Griiner, Rev. Mod. Phys. 66, 1 (1994) and references 
therein. 

69 T. Sasaki and N. Toyota, Synth. Metals, 70, 849 (1995). 

70 F.L. Pratt, T. Sasaki, N. Toyota and K. Nagamine, Phys. 
Rev. Lett. 74, 3892 (1995). 

71 P. Wzietek, F. Creuzet, C. Bourbonnais, D. Jerome, K. 
Bechgaard and P. Batail, J. Physique I 3, 171 (1993). 

72 K. Miyagawa, A. Kawamoto, Y. Nakazawa and K. Kan- 
oda, Phys. Rev. Lett. 75, 1174 (1995). 

73 K. Miyagawa, A. Kawamoto and K. Kanoda, Phys. Rev. 
B 56, R8487 (1997). 

74 N. Harrison, Phys. Rev. Lett. 83, 1395 (1999). 

75 N. Biskup, J. A.A.J. Perenboom, J.S. Brooks and J.S. 
Quails, Solid St. Commun. 107, 503 (1998). 

76 P.L. Kuhns, J.S. Brooks, T. Caldwell, W.G. Moulton, A.P. 
Reyes, N. Biskup, A.M. Kini, J.A. Schlueter, H.H. Wang, 
U. Geiser and J.M. Williams, Solid St. Commun. 109, 637 
(1999). 

77 H. Mori, S. Tanaka, M. Oshima, G. Saito, T. Mori, Y. 



Maruyama, and H. Inokuchi, Bull. Chem. Soc. Jpn. 63, 
2183 (1990). 
L. Ducasse and A. 

(1994). 

H. Seo, |cond-mat/9911329 . 
H. Mayaffre, P. Wzietek, D. 



Fritsch, Solid St. Commun. 91, 201 



Lenoir and P. 



Jerome, C. 

Batail, Phys. Rev. Lett. 75, 4122 (1995). 
A. Kawamoto, K. Miyagawa, Y. Nakazawa and K. Kan- 
oda, Phys. Rev. Lett. 74, 3455 (1995). 
Y. Nakazawa and K. Kanoda, Phys. Rev. B 60, 4263 
(1999). 

T. Komatsu, N. Matsukawa, T. Inoue, and G. Saito, J. 
Phys. Soc. Jpn. 65, 1340 (1996). 

H. Kobayashi, T. Udagawa, H. Tomita, K. Bun, T. Naito 
and A. Kobayashi, Chem. Lett., 1559 (1993). 
H. Tanaka, A. Kobayashi, T. Saito, K. Kawano, T. Naito 
and H. Kobayashi, Adv. Mater. 8, 812 (1996). 
L.K. Montgomery, T. Burgin, J.C. Huffman, J. Ren, and 
M.-H. Whangbo, Physics C 219, 490 (1994). 
H. Kobayashi, H. Akutsu, E. Arai, H. Tanaka and A. 
Kobayashi, Phys. Rev. B 56, R8526 (1997). 
H. Seo and H. Fukuyama, J. Phys. Soc. Jpn. 66, 3352 

(1997) . 

N. Katoh and M. Imada, J. Phys. Soc. Jpn. 63, 4529 
(1994); ibid, 64, 1437 (1995). 

A. Schwartz, M. Dressel, G. Griiner, V. Vescoli, L. De- 
giorgi and T. Giamarchi, Phys. Rev. B 58, 1261, (1998). 
V. Vescoli, L. Degiorgi, W. Henderson, G. Griiner, K.P. 
Starkey and L.K. Montgomery, Science 281 1181 (1998). 
S. Mazumdar, in Interacting Electrons in Reduced Dimen- 
sions, eds. D. Baeriswyl and D.K. Campbell, NATO ASI 
Series B, Vol. 213 (Plenum, NY, 1989), pp. 315 - 329. 
R.L. Greene in Organic Superconductivity, eds. V. Kresin 
and W.A. Little, Plenum, N.Y. (1990), pp. 7 - 13. 
Y.J. Uemura, L.P. Le, G.M. Luke, B.J. Sternlieb, J.H. 
Brewer, T.M. Riseman, G. Saito and H. Yamochi, in ref- 
erence 89, pp. 23 - 29. 

B. Brandow, Phys. Rep. 296, 1 (1998). 
H. Kondo and T. Moriya, J. Phys. Condens. Matter 11, 
L363 (1999). 

H. Kino and H. Kontani, J. Phys. Soc. Jpn., 67, 3691 

(1998) . 

J. Schmalian, Phys. Rev. Lett. 81, 4232 (1998). 
M. Vojta and E. Dagotto, Phys. Rev. B 59, R713 (1999). 
K. Kuroki and H. Aoki, Phys. Rev. B 60, 3060 (1999). 
S. Zhang, J. Carlson and J.E. Gubernatis, Phys. Rev. 
Lett. 78, 4486 (1997). 

M. Guerrero, G. Ortiz and J.E. Gubernatis, Phys. Rev. B 
59, 1706 (1999). 

C. T. Shih, Y.C. Chen, H.Q. Lin and T.K. Lee, Phys. Rev. 
Lett. 81, 1294 (1998). 

E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

D. J. Scalapino, Phys. Rep. 250, 330 (1995). See also, J. 
Low Temp. Phys. 117, 179 (1999). 
G.V. Chester, Phys. Rev. A 2, 256 (1970). 
A.F. Andreev in Progress in Low Temperature Physics, 
edited by D.G. Brewer (North Holland, Amsterdam, 
1982), Vol. VIII. 

108 D. Nelson and M.E. Fisher, Phys. Rev. Lett. 32, 1350 
(1974). 



104 
105 



106 
107 



28 



109 A.J. Leggett, Phys. Rev. Lett. 25, 1543 (1970). 

110 A.S. Alexandres and J. Ranninger, Phys. Rev. B 23, 1796 
(1981); S. Robaskiewicz, R. Micnas and KA. Chao, Phys. 
Rev. B 23, 1447 (1981). Note that these two papers deal 
with the negative U Hubbard model, which leads to super- 
conductivity via a Bose condensation ol small bipolarons. 



As noted in Ref. Ill, this mechanism is unlikely to ap- 
ply to the high temperature superconductors. We stress 
that the effective "negative U, positive V" model we en- 
visage here is quite different from the standard "negative 
U" model in that any bipolarons are expected to be quite 
extended and thus substantially more mobile. See Refs. 



112 for more details. Clearly this issue requires further 



study. 

B.K. Chakraverty, J. Ranninger and D. Feinberg, Phys. 
Rev. Lett. 81, 433 (1998). 

' S. Aubry, p. 217 in Polarons and Bipolarons in High-T c 
Superconductors and Related Materials, edited by E. K. 
H. Salje, A. S. Alexandrov, and W. Y. Liang (Cambridge 
University Press, Cambridge, 1995); L. Proville and S. 
Aubry, Eur. Phys. J. B 11, 41 (1999); L. Proville and S. 
Aubry, Eur. Phys. J. B (to appear). 

; M. Imada, J. Phys. Soc. Jpn. 60, 1877 (1991); Phys. Rev. 
B 48 550 (1993). Note that in the effective spin-Peierls 
model related to our incommensurate BCDW state, we 
do not expect nearest neighbor occupancy of defect pairs, 
as the weak confinement effects arising from the distorted 
lattice will allow spatially extended paired states. 



V.J. Emery ; E. Fradkin and S. Kivelson, cond 



mat/0001077. 



iis M.J. Rice and E.J. Mele, Phys. Rev. B 25, 1339 (1982). 

116 S.C. Zhang, S. Kivelson and A.S. Goldhaber, Phys. Rev. 
Lett. 58, 2134 (1987). 

117 R.T. Clay, S. Mazumdar and D.K. Campbell, unpub- 
lished. 

118 SA. Kivelson, V.J. Emery and H.Q. Lin, Phys. Rev. B 
42, 6523 (1990). 

119 S.M. De Soto, CP. Slichter, A.M. Kini, H.H. Wang, U. 
Geiser and J.M. Williams, Phys. Rev. B 52, 10364 (1995). 

120 Y. Nakazawa and K. Kanoda, Phys. Rev B 55, R8670 
(1997). 

121 S. Belin, K. Behnia and A. Deluzet, Phys. Rev. Lett. 81, 
4728 (1998). 

122 J.M. Schrama, E. Rzepniewski, R.S. Edwards, J. Single- 
ton, A. Ardavan, M. Kurmoo and P. Day, Phys. Rev. Lett. 
83, 3041 (1999). 

123 I.J. Lee, M.J. Naughton, CM. Danner and P.M. Chaikin, 
Phys. Rev. Lett. 78, 3555 (1997). 

124 I.J. Lee, D.S. Chow, W.C Clark, J. Strouse, M.J 



Naughton, P .M. Chaikin and S.E. Brown, cond 



mat /000 1332 



29 



